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ABSTRACT 


Two  separate  non-linear  effects  which  influence  pulsed 
laser  propagation  and  target  interaction  have  been  studied, 
e.g.,  laser  supported  combustion  wave  (LSC)  propagation  and 
laser-induced  gas  breakdown.  A  theoretical  model  describing 
LSC  propagation  and  which  includes  the  non-linear  contribu¬ 
tions  of  the  plasma  transport  and  thermodynamic  coefficients 
throughout  the  wave  profile  has  been  developed.  The  model 
is  one-dimensional  and  the  system  of  equations  is  solved  by 
use  of  the  method  of  quasilinearization.  The  model  utilizes 
the  deflagration  approximation  which  implies  that  pressure 
balance  with  the  atmosphere  exists.  Numerical  results  of  our 
calculations  are  in  good  quantitative  agreement  with  available 
experimental  data.  Further  numerical  calculations  have  been 
performed  by  use  of  a  one-dimensional,  variable  area  hydro- 
dynamic  code.  It  is  found  that  a  precursor  shock  wave  which 
is  driven  by  the  LSC  wave  precedes  the  LSC  wave  and  that  attain¬ 
ment  of  pressure  balance  with  the  atmosphere  does  not  occur  on 
a  time  and  spatial  scale  consistent  with  experimental  observa¬ 
tions.  These  calculations  have  been  performed  including  radial 
mass  loss  and  beam  divergence  effects.  It  is  postulated  that 
two-dimensional  flow  effects  in  the  vicinity  of  the  target, 
which  have  not  been  specifically  included  in  the  calculation 
may  be  important.  Inclusion  of  such  effects  on  an  ad  noc 
basis  results  in  good  correlation  with  experimental  data. 

Experiments  have  been  performed  to  determine  the  effect 
of  particulate  matter  on  laser-induced  air  breakdown  for  Nd 
laser  wavelengths  (1.06y)  and  for  pulse  widths  of  30  nsec. 

We  have  measured  breakdown  thresholds  for  air  in  the  presence 
of  Si02  and  carbon  fibers  of  9y,  6.6y  and  68y  diameters.  It 
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1.  INTRODUCTION  AND  SUMMARY 

The  work  described  in  this  report  was  performed  at 
Systems,  Science  and  Software  (S3)  during  the  period 
21  December  1973  through  21  June  1974  on  Contract  No. 

DAAH01-74-C-0386 .  The  objective  of  the  work  is  twofold; 

(1)  to  develop  a  theoretical  model  including  a  description 
of  the  various  non-linear  effects  for  laser  supported  com¬ 
bustion  wave  (LSC )  propagation,  and  (2)  to  perform  measure¬ 
ments  at  1.06 vi  to  determine  the  effect  of  particulate  matter 
on  laser-induced  air  breakdown.  The  latter  work  was  performed 
by  S3  personnel  at  the  Battelle  Columbus  Laboratories  under 
a  subcontract  from  S3. 

Previous  work^1  has  shown  the  existence  of  LSC 

waves  initiated  at  target  surfaces  for  laser  fluxes  as  low  as 
4  2 

4><io  watts/cm  for  CC>2  laser  radiation.  The  existence  of 

such  wave  phenomena  has  important  implications  regarding  the 

delivery  of  laser  radiation  directly  to  the  surface  since  the 

r  4 1 

waves  are  highly  absorptive.  Raizer  argued  that  the  propa¬ 
gation  and  structure  of  such  waves  could  be  adequately 
described  by  a  one-dimensional  model  which  included  propaga¬ 
tion  of  the  wave  away  from  the  surface  by  thermal  conduction. 

He  utilized  a  linearized  analysis  of  the  gasdynamics  conser¬ 
vation  equations  to  obtain  estimates  for  the  wave  velocity 

as  a  function  of  the  incident  laser  flux.  The  work  at  Mathe- 

r  3  ] 

matical  Sciences  Northwest  showed  that  the  LSC  wave  velo¬ 
city  was  much  greater  (about  an  order  of  magnitude)  than  could 
be  predicted  from  Raizer 's  theory,  and  that  the  velocity  did 
not  scale  with  flux  according  to  Raizer* s  prescription.  Sub¬ 
sequently,  we  pointed  out  that  in  addition  to  thermal  conduc¬ 
tion,  the  transport  of  thermal  radiation  in  the  vacuum  ultra¬ 
violet  region  was  the  dominant  vrave  transport  mechanism.  ^ 4  *  ^ 
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Indeed,  the  order  of  magnitude  discrepancy  with  experiment 
was  removed.  However,  our  analysis  was  also  based  on  a 
linearized  model.  The  present  work,  which  is  an  extension 
of  of  our  previous  work,  removes  the  linearized  assumptions 
and  includes  the  non-linear  variation  of  transport  and  thermo¬ 
dynamic  coefficients  in  the  analysis.  Two  separate  analyses 
are  presented  in  Chapters  2  and  3.  In  Chapter  2,  we  utilize 
the  method  of  quasilinearization  to  solve  the  LSC  propagation 
problem.  it  is  shown  that  the  non-linear  features  of  the 
problem  are  profound  and  that  by  taking  them  into  account, 
good  agreement  is  found  with  the  MSNW[3]  experimental  results. 

However,  this  work  utilizes  the  deflagration  approxima¬ 
tion,  i.e.,  the  pressure  ahead  of  and  through  the  wave  is 
approximated  as  a  constant.  To  verify  this  assumption  and  to 
study  the  importance  of  potential  2-D  effects,  we  describe  in 
Chapter  3  some  numerical  calculations  based  on  a  one-dimensional 
variable-area  hydrodynamics  code.  These  calculations  indicate 
that  a  precursor  shock  wave  exists  ahead  of  the  LSC  wave,  thus 
jeopardizing  the  validity  of  the  deflagration  approximation. 

It  is  shown  that  radial  expansion  of  the  laser-heated  plasma 
beyond  the  confines  of  the  laser  beam  is  not  a  significant 
loss  mechanism.  However,  we  do  find  that  the  presence  of  the 
target  results  in  higher  wave  velocities  than  are  predicted 
by  previous  analyses.  These  velocities  are  higher  than  those 
observed  experimentally [3 ]  at  similar  distances  from  the  tar¬ 
get.  We  conclude  that  mass  loss  from  the  plasma  in  the  vicinity 
of  the  target  must  be  responsible  for  this  more  rapid  attainment 
of  pressure  balance  with  the  atmosphere.  Upon  presuming  on  an 
ad  hoc  basis  the  presence  of  such  a  loss  mechanism,  our  calcu¬ 
lations  are  in  good  agreement  with  the  MSNW  data  as  well  as 
with  those  results  described  in  Chapter  2. 
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Finally,  in  Chapter  4,  we  discuss  measurements  per¬ 
formed  by  us  at  1.06y  to  determine  the  wavelength  and  pulse 
length  scaling  which  exists  between  10.6m  and  1.06y  and  for 
pulse  length  variations  between  300  nsec  and  30  nsec.  We  per¬ 
formed  measurements  on  Sit^  and  carbon  fibers  in  an  air  environ¬ 
ment  and  find  significant  reductions  in  the  laser-induced  break¬ 
down  threshold  from  the  clean  air  value.  For  example,  we  find 

for  a  68yi  carbon  fiber  a  reduction  in  breakdown  threshold  to 
10  2 

2.3xio  watts/cm  .  A  full  order  of  magnitude  less  than  the 
clean  air  value. 

It  is  shown  that  the  breakdown  phenomena  at  both  10.6m 
and  1.06m  are  similar  and  that  the  breakdown  thresholds  (for 
a  particular  pulse  length)  scale  as  the  inverse  square  of  the 
wavelength.  This  conclusion  is  consistent  with  arguments 
based  on  inverse  bremsstrahlung  as  the  dominant  absorption 
process. 
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2 .  NONLINEAR  THEORY  OF  LASER  SUPPORTED  COMBUSTION  WAVES 

2 . 1  INTRODUCTION 

It  has  been  observed  experimentally^1'2'3^  that  a  sub¬ 
sonic  wave  can  propagate  and  is  supported  by  the  addition  of 
energy,  e.g.,  laser  radiation.  The  earlier  theoretical  work 
on  laser  supported  combustion  wave  (LSC)  propagation  was 
based  on  a  one-dimensional  linear  theory  including  thermal 
conduction  as  the  only  energy  transport  mechanism. ^  Later, 
it  was  pointed  out  by  Boni  and  Su  that  thermal  radiation  in 
the  vacuum  ultraviolet  (VUV)  spectrum  is  the  dominant  energy 
transfer  mechanism.  Depending  upon  the  absorption  coeffi¬ 

cient  in  the  VUV  band,  we  have  demonstrated  that  the  wave 
velocity  can  vary  linearly  (when  the  vacuum  ultraviolet  pre¬ 
cursor  is  not  optically  thick)  with  the  incident  flux  or  as 
a  square  root  function  of  the  incident  laser  flux  (when  the 
precursor  can  be  considered  as  optically  thick) . ^  However, 
all  the  above  results  are  based  on  analyses  which  result 
from  a  linearization  of  the  governing  equations.  In  order  to 
obtain  a  better  understanding  of  the  interaction  between  the 
laser  radiation  and  gasdynamic  processes  involved  in  LSC  pro¬ 
pagation,  a  detailed  treatment  of  the  problem,  which  includes 
a  description  of  the  various  non-linear  effects  is  required. 

In  this  report,  we  first  set  up  a  nonlinear  mathematical  model 
for  a  one-dimensional  steady  state  LSC  wave  propagating  into 
air  with  real  gasdynamic  properties.  Then  we  proceed  to  solve 
the  problem  by  use  of  the  method  of  quasilinearization, 
which  has  been  used  quite  extensively  in  boundary  layer 
theory . 1(^  Since  the  equation  is  nonlinear,  there  are 
certain  difficulties  associated  with  the  process  of  obtaining 
a  convergent  solution.  We  discuss  these  difficulties  and 
finally  present  some  solutions  and  compare  with  available 
experimental  data. 
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2  •  2  FORMULATION  OF  TIIE  PROBLEM 

Assume  a  parallel,  axially-symmetric  laser  beam  with  an 
intensity  lower  than  the  breakdown  flux  (<_  107  watts/cm2)  inci¬ 
dent  onto  a  partially  absorbing  gas.  Since  the  laser  radiation 
cannot  support  a  laser  supported  detonation  (LSD)  wave,  any 
initial  shock  waves  present  due  to  target  interaction  have  pro¬ 
pagated  far  beyond  the  light-absorbing  plasma.  Therefore,  we 
assume  that  the  pressure  in  the  absorption  wave  has  equilibrated 
with  the  ambient  and  a  constant  pressure  process  is  obtained. 

although  it  has  been  suggested  that  two-dimensional  ef¬ 
fects  are  important,  in  this  study  we  assume  a  one-dimension 
planar  geometry  in  order  to  simplify  the  problem.  However, 
the  heat  diffusion  through  the  radial  direction  will  be  approxi¬ 
mately  included  as  Raizer  suggested. in  Section  3,  the  two- 
dimensional  (2-D)  effects  are  studied  by  using  a  quasi  one¬ 
dimensional  approach  and  discussed. 

We  assume  for  the  present  that  a  continuous  laser 
source  is  used  such  that  the  LSC  can  propagate  at  a  steady 
speed,  vq,  into  the  ambient  gas  with  density  pq  and  tempera¬ 
ture  Tq.  Accordingly,  coordinates  are  transformed  into 
the  wave  frame  (Figure  1) .  in  this  coordinate  system,  the 
continuity  equations  can  be  simply  integrated  to  give: 


where  p  and  u  are  the  density  and  local  velocity  at  any  posi¬ 
tion  in  the  wave.  The  energy  equation  can  also  be  written  as 


(2) 
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where  h  is  the  specific  enthalpy  and  X  is  the  thermal  conduc- 

c 

livity.  In  the  above  equation,  radiative  transport  is  in¬ 
cluded  through  the  terras  which  is  the  integrated  flux,  and 
4>  which  is  the  incident  laser  flux.  The  term  j  represents 
the  volumetric  energy  loss  through  the  lateral  surface  of  the 
light  channel.  Assuming  the  plasma  in  the  light  channel  to 
behave  like  a  heated  cylinder  with  strongly  cooled  walls,  one 
can  show  that  the  two  dimensional  heat  loss  by  conduction 
through  the  channel  walls  can  be  approximated  as 


j 


c 


f 


(3) 


where  A  is  a  constant  of  order  of  unity. 

The  volumetric  laser  heating  term  obeys  the  following 
transport  equation; 


dcj) 

dx 


=3 


k,  $ 

Li 


with  the  boundary  condition  4>  =  4>  at  x  =  and  kT  is  the 

-1  L 

linear  absorption  coefficient  (cm  ) . 

The  are  three  boundary  conditions  which  the  system  of 
equations  (1)  and  (2)  must  satisfy.  At  the  far  upstream 
boundary,  the  gas  is  cold  at  a  temperature  Tq,  and  also  the 
thermal  flux  vanishes,  i.e., 


T 

(i 


dT 

dx  - 


0 


at  x  = 


(4) 


(5) 
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On  the  other  hand,  far  downstream,  the  energy  gain  through 
absorption  of  the  laser  energy  must  be  balanced  by  the  energy 
loss  through  radiation  and  conduction;  therefore 

dT  . 

37  “  0  at  x  =  »  (6) 

Hence  the  problem  becomes  an  eigenvalue  problem.  A  unique 
solution  exists  only  at  a  discrete  value  of  the  velocity,  v  , 
which  is  the  eigenvalue.  This  velocity  is  the  one  at  which° 
the  wave  propagates. 


2.3  RADIATIVE  AND  GASDYNAMIC  PROPERTIES  OF  THE  PLASI1A 

In  a  high-temperature  plasma,  the  absorption  coefficient 
of  the  laser  light  with  frequency  v  is  dominated  by  the  free- 
free  bromsstrahlung  radiation  of  electrons  in  both  the  ionic 
and  atomic  fields.  The  absorption  coefficient  for  ionic 
bi emsstrahlung  may  be  calculated  by  the  Kramer's  formula, 
corrected  for  stimulated  emission  given  below, 


Here  G  is  the  Gaunt  factor 


Assuming  for  the  present  that  the  state  of  ionization  may  be 
described  by  the  Saha  equation,  the  number  density  of 

electrons  (Ne)  or  ions  (H+)  can  then  be  calculated  according 
to 
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where  J  is  the  total  number  density  (IJ  =  P/kT)  ,  where  P  is  the 
total  pressure  due  to  atoms,  ions,  and  electrons.  The  ratio 

g+/ga  is  approximately  1.9  for  air,  whereas  the  averaged 
ionization  potential  is  given  as  14.4  eV. 

As  we  have  demonstrated  before, [5]  thermal  radiation 
processes  are  very  important  for  the  gas  in  the  temperature 
range  from  1  eV  to  2  eV.  Since  radiative  properties  are  fre¬ 
quency  dependent,  it  is  common  practice  to  use  a  step  model 
(or  a  band  model)  in  most  of  the  high  temperature  gasdynamic 
calculations.  Such  a  model  consists  of  a  different  constant 
mean  absorption  coefficient  in  various  frequency  ranges.  In 
this  report,  for  simplicity,  we  adopt  the  two-step  absorption 
coefficient  model  (Figure  2)  to  investigate  the  nonlinear 
effects  of  LSC  wave  propagation.  Accordingly,  the  integrated 
radiative  flux  can  be  separated  into  two  parts,  F  =  F  +F  , 
where  covers  the  visible  and  infrared  portions  of  the2 
spectrum  ,  and  F 2  covers  the  vacuum  ultraviolet  (VUV)  portion 

of  the  spectrum.  The  linear  absorption  coefficients,  k  and 

k2'  by  usin9  Allen's  data [ 12 ]  can  be  simply  correlated  as[13] 
(at  p  =  1  atm) 


kl  =  °-00255  fa) 


4.34 


cm 


-1 


1.776 


fa) 


-0.71 


cm 


-1 


1130  A  <  X  <  co 


0  <  X  <  1130  A 


(9) 


(10) 


It  can  be  seen  that  for  a  wave  thickness  of  a  few  centimeters 
the  long  wavelength  radiation  can  be  considered  as  transparent, 
whereas  for  the  VUV  portion,  the  gas  is  not  strictly  optically 
thick.  However,  it  can  be  shown  by  a  simple  calculation  based 
on  the  exact  integral  expressions,  that  the  flux  calculated  by 
the  optically  thick  band  approximation  overestimates  the  radia¬ 
tive  flux  by  at  most  a  factor  of  two  in  the  temperature  range 


9 


Schematic  diagram  of  a  two  band  absorption  coeffi 
cient  model. 


R-2344 


oL  interest.  it  is  our  purpose  in  the  present  report  to 
•love lop  a  simple  numerical  scheme  in  solving  the  LSC  wave 
problem.  Therefore,  the  simple  optically  thick  band  approxi¬ 
mation  is  adopted  in  this  report.  Hence  the  radiative  flux 
divergence  in  Equation  (2)  can  be  expressed  as[5] 


dF 


ur  ,  , 

— £  =  _  1_  U  dT  \  . 

dx  dx  \Ar  ZZ)  + 

where  Ar  is  the  radiative  conductivity  defined  as: 


(11) 


X  =  4 _  ^2 

r  3k2  dT 


aT3  “u 
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(12) 


Also  U;L  -.hV-j/kT  and  corresponds  to  the  wavelength 

X  =  113  0  A.  The  —  Jt  _  .  .. 


A  1130  A.  The  term  jr  represents  the  transparent  cooling 


part.  Since  95  percent  of  the  blackbody  energy  falls  within 
the  visible  and  IR  spectrum  for  the  temperature  range  of 
interest,  we  simply  approximate  the  cooling  term  as 


jr  =  4irk 


if  D 


BJT)  -  BV(TJ  dv 


-  4ak 


1 (t4-Tc ) 


(13) 


For  the  thermodynamic  properties  of  high  temperature 
air,  it  has  been  demonstrated  by  Su[14J  that  a  simple  power 
law  correlation  can  retain  both  reality  and  simplicity.  Under 

a  constant  pressure  process  (p  =  1  am),  the  enthalpy  can  be 
correlated  as 


h  =  5.036  x  10 
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f  151 

from  Predvoditelev* s  data  (Figure  3).  The  thermal  conduc¬ 
tivity  of  high  temperature  air  has  been  calculated  by 

r  1 6 1 

Hanson.  As  one  can  see  (Figure  4),  the  data  are  highly 

oscillatory  because  of  the  changes  of  chemical  composition 
with  temperature.  However,  because  of  the  range  of  interest 
(1  eV  >  T  >  2  eV)  and  simplicity,  calculated  X^  (according 
to  Equation  12)  and  X^  are  combined  and  correlated  to  two 
single  exponential  functions  as  shown  in  Figure  4, 

„  1vln8  6 . 4  (T/104-l  .75)  , 

X (T)  =  X  +  X  =  2.2x10  e  '  erg/sec  cm  °K) 

J L  O 

T  <  17500  °K 

_  _v1n8  3 (T/104-1.75)  , 

=  2.2x10  e  erg/sec  cm  °K 

T  >  17500  °IC  .  (14) 


On  the  other  hand,  in  the  conductive  energy  loss  term  j 
Equation  (3),  X c  is  simply  correlated  by  the  following  for¬ 
mula  as  in  Reference  17  (Figure  5) : 


2.4  THEORY  OF  QUA3ILIHEARIZATI0N 

Because  of  the  inherent  nonlinearity  of  the  LSC  wave 
problem,  as  described  above  one  has  to  resort  to  numerical 
integration  of  the  governing  differential  equations.  The 
principal  difficulty  in  numerical  analysis  of  this  problem 
is  the  split,  two-point  nature  of  the  boundary  conditions. 
Since  initial  conditions  are  required  to  start  the  numerical 
integration,  a  guess  based  on  physical  intuition  of  the  ini¬ 
tial  condition  has  to  be  made.  In  general,  however,  after 
the  numerical  integration,  the  boundary  conditions  on  the 
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end  point  are  not  satisfied.  In  order  to  satisfy  the  second 
boundary  condition,  one  must  iterate  until  convergence,  within 
some  predetermined  accuracy  criterion,  is  obtained.  A  common 
method  is  to  obtain  additional  solutions  with  different  ini¬ 
tial  guesses,  followed  by  inverse  interpolation  or  extrapola¬ 
tion  on  the  boundary  value  attained  at  the  end  point.  Another 
method,  called  successive  substitution,  involves  solving  for 
the  highest  derivative  and  substituting  the  results  successive¬ 
ly  into  the  right  hand  side  starting  with  an  approximation 
that  satisfies  the  boundary  conditions.  Both  methods  require 
a  reasonable  approximation  to  the  boundary  condition  in  the 
first  method  and  functions  in  the  second  method.  Without  a 
sufficiently  good  initial  guess,  the  processes  may  diverge  or 
converge  very  slowly.  Such  a  technique  is  not,  in  general, 
desirable . 

Quasilinearization  may  be  viewed  as  an  extension  of 
Newton' s  method  for  the  solution  of  algebraic  equations  to  the 
solution  of  differential  equations . [ 7]  This  method  has  been 
widely  used  in  solving  boundary  layer  problems ' 9 ' 10 ^  and 
has  been  demonstrated  to  be  a  powerful  method  for  handling 
split  point  boundary  conditions. 

Consider  a  system  of  nonlinear  first  order  ordinary 
differential  equations: 

dyi 

=  fi<Yl.y2"-yn)  i  =  1,2*  *  *n  ,  (15) 

where  is  the  ith  dependent  variable,  and  x  is  the  indepen¬ 
dent  variable,  where  some  of  the  boundary  conditions  speci¬ 
fied  at  x  =  0, 

y^O)  =  Bi  i  =  1,2  •••*  ,  (16) 

and  some  are  specified  at  the  other  boundary,  x  =  a, 

y^(a)  =  Bi  i  =  Jl+1 ,  •  •  *n  .  (17) 
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"ho  (k+1)  approximation  to  the  solution,  y.(k+1)  ,  can 
be  obtained  by  expanding  the  right-hand  side  of  the  equation 

about  the  kth  approximation  and  retaining  only  the  linear 
terms  : 


dy 


k+1 


dx 


=  qiy*) 


t  y--(yk)  ry^ 


This  equation  is  linear  in  yk+1;  and  the  system  of  equations 
can  be  solved  by  any  standard  numerical  integration  routine. 
Also,  since  it  is  a  linear  equation,  the  principal  of  super¬ 
position  can  be  applied. 


A  particular  solution,  P . (x) ,  may  be  obtained  by  inte¬ 
grating  the  entire  equation,  starting  with  the  initial  condi- 
tions : 


Yj.  (0)  =  Pi<0)  =0  i  =  l/2/- -n 


(19) 


Then  one  homogeneous  solution  is  obtained  for  each  of  the  n 
boundary  conditions  with  the  initial  conditions: 


yi(0)  -  H..(0)  =  6.. 


i/j  —  1,2, 


•  n 


(20) 


where  6±.  is  the  ICronecker  delta,  6..  =  i  when  i  =  j  and 

6ij  =  °  When  1  ?  >  The  requirement  that  the  sums  of  the 
particular  and  homogeneous  solutions  satisfy  the  boundary 

conditions  determines  the  combination  constant  C  • 

j  ‘ 


n 

g  CjVV  -  Bi-pi(bi; 


i  1/2 , 


•n 


(21) 


v/iiere  b  •  =0  for  i  =  1  ?  •«.p  k 

l  ^  and  bi  =  a  for  i  =  £+!••. n. 
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-he  new  approximation,  y  ,  is  then  calculated  at  each 
point  by  summing  the  stored  values  of  P  and  II 


n 


y^+1(x)  =  p.  (x)  +  2  C.  II. 
1  1  j=i  3  i 


;(x) 


(22) 


or  by  obtaining  the  initial  conditions  in  this  way  and  inte¬ 
grating  forward. 

notice  that  if  the  particular  solution,  P,  is  required 
to  satisfy  the  i  known  initial  conditions,  then  there  are 
only  (n-£ )  remaining  boundary  conditions  to  satisfy,  or  there 
are  only  (n— £ )  homogeneous  equations  to  solve,  which  dras¬ 
tically  reduces  the  integration  time. 

Combining  Equations  (1),  (2)  and  (11)  gives  a  second 

order  nonlinear  equation  governing  the  propagation  and  struc- 
tue  of  the  LSC  wave. 


d 

as 


P 


(22a) 


I 

Here  X^,  X^,  h,  kL<f ,  jg  and  jr  are  all  complicated,  non-linear 
functions  of  T.  Defining  X  =  X C+Ar  and  j  =  jc+jr,  it  is  shown 
in  Appendix  A  that  Equation  (22a)  can  be  rewritten  as 


d2Tk+1 

dx2 


k  dT  f 
P  dxi" 


+1 


_ rkrj?k+l  +  sk  _ 


(23) 


for  -»  <  x  <  +  oo,  vnere 


Sk  =  +  ,  V 


Hk 


(24) 


Details  of  the  method  of  solution  and  it’s  associated  diffi¬ 
culties  appear  in  App<ndix  A.  The  interested  reader  should 
refer  to  that  section  of  the  report  for  further  discussion 
of  the  implementation  of  the  method  of  quasilinearization  and 
its  intricacies. 
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2  •  RESULTS  AND  DISCUSSION 

Before  we  discuss  the  results  of  nonlinear  calcula¬ 
tions,  let  us  briefly  describe  some  of  the  linear  results 
from  References  6  and  17  which  will  be  used  for  comparison 
For  a  linear  model.  Equation  (2)  can  be  simplified  as 


d2?  _  PoVoCp  dT  G  „ 
dx2  X  dx  "  X  1 


=  0  for  T  <  T. 

in 


x  <  0 


with  boundary  conditions 


(25) 


and 


at  x  =s  ±oo 


T  -  0  at  x  =  -oo 


All  the  coefficients  in  Equation  (25)  are  constants.  Compar¬ 
ing  Equation  (25)  with  Equation  (22a),  one  recognizes  that 
CpT  is  equivalent  to  k(T)  and  GT  is  equivalent  to 

+  3r(T)«  Since  Equation  (25)  is  a  linear  equation, 
it  can  be  solved  analytically  to  give 


(26) 

where  is  the  ignition  temperature  at  x  =  0  (usually 

Tin  =  12 ' 000  °K>*  The  temperature  profile  is  given  by  the 
equations 


T  =  T.  e 
m 


“a1 I x I 


x  <  0 


a0x  k_4>  /A 

T  =  T._  e  2  +  -L-  -°  ,  /  a2x  \ 

ala2  \e  “  1J  x  >  0  , 


(27) 
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Again,  to  present  some  of  the  nonlinear  calculations, 
we  perform  calculations  for  the  experiment  of  Bunkin,  et  al.^ 
"hey  indicate  that  for  an  experiment  utilizing  Nd  laser- 
irradiated  air,  the  temperature  (T^)  in  the  downstream  portion 
of  the  wave  is  estimated  around  18,000  °K  for  p  =  1  atm.  The 
laser  power  was  1,100  kw  incident  onto  a  0.25  cm  spot  which 
yields  a  flux  of  <J>o  —  2.25x10  watts/cm3.  The  measured  laser 
absorption  coefficient  is  about  7xio“3  cm-1.  Therefore  in 
the  downstream  portion  of  the  wave,  the  absorbed  radiative 
energy  can  be  approximated  as  kL<j>  «  1.57xl05  watts/cm3,  where 
in  the  upstream,  practically  no  energy  is  absorbed  (or  kL  -  0) . 
Using  averaged  values  of  X,  and  G,  the  energy  losses 
through  radiation  and  conduction  account  for  21.6%  of  the 
total  energy  absorbed.  With  these  losses  the  velocity  pre¬ 
dicted  by  the  linear  theory  is  about  9.36  m/sec.  [ 61  The  tem¬ 
perature  profile  is  then  calculated  from  Eq.  (27)  and  plotted 
in  Figure  6.  The  wave  front  thickness  is  greater  than  10  era 
and  the  temperature  reaches  54,000  °K  (this  can  be  seen  simply 
by  letting  dT/dx  =  0,  and  then  Tro  =  kL<j»0/G)  .  Now  let  us 
examine  the  nonlinear  effects  on  the  LSC  wave  by  gradually 
replacing  the  linear  coefficients  in  Eq.  (25)  with  the  corres¬ 
ponding  nonlinear  terms.  By  introducing  the  real  equation  of 
state  (curve  B  in  Figure  6),  the  wave  front  spreads  even  more 
(final  temperature  reaches  54,000  °K  at  a  distance  of  about 
35cm  downstream)  ,  due  to  the  fact  that  most  of  the  energy 
absorbed  in  the  wave  front  is  used  to  ionize  the  gas  first 
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instead  of  increasing  the  temperature  linearly  as  in  the  ideal 
gas  law.  next,  the  radiative  energy  loss  in  the  linear  equa¬ 
tion  is  replaced  by  the  real  loss?  Eq.  (13).  This  has  a  tre¬ 
mendous  influence  on  the  temperature  profile  (curve  G  in 
Figure  6)  since  the  real  radiative  loss  term  is  proportionate 
to  T  rather  than  T.  This  effect  reduces  the  downstream  tem¬ 
perature  from  54,000  °K  to  about  26,000  °K,  and  at  the  same 
time  the  wave  front  thickness  is  reduced  to  4  cm.  The  tem¬ 
perature  profile  steepens  even  more  as  the  nonlinear  diffusivity 
[X  =  AC(T)  +  Ar(T)]  is  introduced  (curve  D  in  Figure  6).  It 
can  be  seen  that  a  definite  wave  front  is  formed  within  about 
2  cm.  Finally,  with  the  nonlinear  conductive  loss  included 
(curve  E  in  Figure  6),  the  downstream  temperature  drops  even 
further,  to  a  reasonable  range  T^  «  19,400  °IC.  From  the  above 
discussion,  it  is  clear  that  the  nonlinear  radiative  and  con¬ 
ductive  diffusivity  establish  the  shape  of  the  wave  front 
whereas  the  radiative  and  conductive  loss  mechanisms  deter¬ 
mine  the  final  temperature  profile  behind  the  wave.  The 
wave  speed  predicted  by  the  above  cases  varies  from  4.0  m/sec 
(case  E)  to  13.2  m/sec  (case  D) .  We  will  discuss  more  about 
the  wave  velocity  in  the  next  few  paragraphs. 

Figures  7  and  8  compare  the  energy  balance  for  the 
following  three  cases  a)  kL4>o  =  2.5xi05  watts/cm3  and 
R  =  0.125  cm  [where  R  is  the  radius  of  the  spot,  Eq.  (14)], 

b)  ^o  =  1*57x105  watts/cm3  and  R  =  0.125  cm  and 

c)  Vo  =  1-57x105  watts/cm3  and  R  =  0.25  cm.  We  plot  the 
temperature  profiles  of  the  above  three  cases  in  Figure  7. 

The  wave  fronts  are  almost  exactly  the  same  whereas  the  down¬ 
stream  temperatures  in  the  three  cases  are  different.  With 
high  absorbed  energy,  Case  A  shows  a  higher  temperature  than 
Case  B.  However,  with  the  different  spot  size,  the  tempera¬ 
ture  of  the  LSC  wave  with  lower  energy  addition  (Curve  C)  can 
be  higher  than  the  temperature  profile  with  larger  energy 
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Figure  8  Energy  balance  for  three  different  cases.  (A) 
2.5  x  10  watts/cm3  and  R  =  0.125  cm.  (B)  k2c J> 
1.57  x  10 5  watts/cm3  and  R  =  0.125  cm  and  (C) 
k2<f>  -  1.57  x  10 5  watts/cm3  and  R  =  0.25  cm. 
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addition  (Curve  A).  It  is  not  surprising,  since  the  conduc¬ 
tive  energy  is  inversely  proportional  to  the  spot  size,  i.e., 
jr  a  1/R  .  Therefore  in  Case  C  the  conductive  loss  is  only 
about  one  fourth  of  the  loss  in  Case  A.  Figure  3  shows  the 
comparison  of  absorbed  energy,  radiative  loss  term  and  con¬ 
ductive  loss  term  for  all  three  cases.  For  Cases  A  and  B 
the  conductive  losses  are  comparable  while  the  radiative 
losses  show  a  tremendous  difference.  This  effect  is  in 
part  due  to  the  smaller  thermal  conductivity  at  higher  tem¬ 
perature  (see  Figure  4)  and  the  strong  temperature  depen¬ 
dence  of  the  radiative  energy  loss.  The  absorbed  energy  in 
Cases  B  and  C  are  the  same,  however,  Case  B  is  dominated  by 
conduction  loss  and  radiative  loss  dominates  for  Case  C.  It 
is  interesting  to  notice  that  the  area  of  -  j  -j  of 

Case  A  and  C  are  about  the  same.  This  is  the°reason  Shy  the 
temperature  profile  are  about  the  same  for  Case  A  and  C 
shown  in  Figure  7,  even  though  their  dominant  loss  mechanism 
is  different.  Since  the  wave  velocity  depends  ppon  the 
amount  of  energy  transported  from  downstream  to  upstream,  we 
expect  that  the  wave  speeds  for  Case  A  and  C  will  be  similar. 
In  fact,  we  obtain  vq  =  10.2  m/sec  for  Case  C  and  v  =9.3 
m/sec  for  Case  A.  Notice  that  the  slope  of  Curve  C  in 
Figure  7  is  slightly  greater  than  the  slope  of  A  therefore 
the  wave  speed  of  C  is  also  slightly  higher  despite  the 
larger  absorbed  energy  in  Case  A.  Curve  C  has  a  much  slower 

velocity;  vq  «  4  m/sec  due  to  smaller  gradient  at  the  wave 
front. 

Figure  9  summarizes  our  calculations  for  the  wave 
velocity  as  a  function  of  absorbed  energy  with  various  beam 
radii.  As  R  decreases,  conductive  energy  loss  through  the 
radial  direction  increases.  When  jr  becomes  the  dominant 
loss  mechanism  the  wave  velocity  decreases  drastically  with 
decreasing  radius.  This  effect  can  be  seen  more  clearly  in 


R-2344 

Figure  10  where  wave  speed  is  plotted  as  a  function  of  radius 
with  various  absorbed  energy  levels  as  a  parameter .  When  the 
spot  size  is  larger  than  one  centimeter  (or  R  >  0.5  cm)  the 
beam  radius  has  little  effect  on  the  wave  propagation  speed. 

In  fact,  for  all  practical  purposes,  the  velocity  can  be  pre¬ 
dicted  without  including  the  conductive  energy  loss  for 
R  >  0.5  cm  in  most  cases  (except  at  very  low  .  The  rea¬ 

son  for  neglecting  j  for  R  >  0.5  cm  is  obviously  shown  in 
Figure  8,  Case  C  where  the  radiative  energy  loss  term  dominates 
the  conductive  energy  loss  even  at  R  =  0.25.  In  the  range 
for  R  <  0.5,  in  general  the  velocity  varies  from  linear  to 
square  term  of  radius  (Figure  10) .  It  is  also  interesting  to 
note  that  for  most  of  the  cases,  the  wave  speed  is  neither  in 

the  linear  regime  (v  a  kT 4  ),  nor  square  root  regime 

o  Jj  o  .  75 

(v  a  /kL4>Q)  »  but  is  proportional  to  tkL<t>0)  * 

For  comparison,  the  prediction  of  wave  speed  from  a 
linear  model  is  also  plotted  in  Figure  9.  It  is  shown  that 
the  wave  speed  has  a  square  root  dependence  on  kL<J>Q  if  kLc|>o 
is  large,  and  a  linear  dependence  close  to  the  threshold  as 
previously  discussed.  Judging  from  the  simplicity  of  the 
linear  theory,  the  prediction  of  wave  speed  has  a  reasonable 
agreement  in  comparison  with  the  nonlinear  theory,  except 
for  the  cases  where  R  is  very  small.  However,  as  we  pointed 
out  before, ^  a  correct  value  of  R  is  not  very  well  defined 
in  a  focused  laser  beam.  Conductive  loss  calculated  in 
Eq.  (3)  is  strictly  applied  only  to  cylindrical-like  plasma. 
Therefore,  for  a  simple  estimate  of  wave  propagation  velocity 
it  is  our  opinion  that  either  a  phenomenological  analysis  or 
a  linear  model  is  sufficient. 

So  far  we  have  not  included  in  our  analysis  the  non¬ 
linear  variation  of  absorption  coefficient,  and  have  treated 
the  cases  where  absorbed  flux  can  be  approximated  as  kL<|>  . 
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It  is  unfortunate  that  we  have  found  that  there  exists  a 
strong  singularity  downstream  if  the  flux  at  any  point  is 
treated  as 


-r 

—  CO 


k (T) dx 


<j>  (x)  =  <j>  e 


As  the  LSC  wave  propagates  the  incident  flux  is  absorbed  by 

the  wave  front,  and  therefore  attenuated  by  the  hot  plasma. 

When  most  of  the  incoming  energy  is  absorbed,  the  temperature 

at  the  far  downstream  is  expected  to  decrease  as  shown  in 

Figure  1.  Hence  there  is  a  maximum  temperature  existing  in 

the  downstream.  At  that  maximum  temperature  the  gradient  of 

the  temperature  profile  is  zero;  dT/dx|  =  0.  Then  the 

*m 

derivative  of  <J>  with  respect  to  T  (where  T  is  treated  as  an 

independent  variable  in  the  method  of  quasilinearization)  is 

singular  at  s  ;  i.e., 
m 


d_ 

dT 


xm 


d  dx 
dx  dT 


m 


Hence  the  method  of  quasilinearization  fails  to  apply  in  the 

region  around  x  . 

m 

Although  the  nonlinear  absorption  effect  has  not  been 
included,  we  feel  that  we  have  already  demonstrated  the  impor¬ 
tance  of  energy  loss  mechanisms  and  the  diffusive  effects  in 
predicting  the  temperature  profile  and  wave  velocity  of  a 

LSC  wave.  We  feel  that  the  parameter  kT  <j>  (or  the  total 

i-i  o 

absorbed  flux)  is  a  useful  parameter  in  interpreting  any 

experimental  data.  For  example  the  data  of  LSC  experiments 

[31 

performed  by  Klosterman,  e_t  al.,  using  a  CO  2  laser  inci¬ 
dent  on  an  A2.  target  in  air  has  been  plotted  on  Figure  9  with 
an  estimated  absorption  coefficient  of  order  unitv.  It  shows 
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very  good  agreeraent  with  the  non-linear  theory,  especially 
the  dependence  on  the  absorbed  flux. 

In  conclusion,  we  have  used  the  method  of  quasilineari¬ 
zation  to  study  nonlinear  effects  of  a  propagating  LSC  wave. 

It  is  found  that  nonlinear  effects  are  very  important  in  pre¬ 
dicting  the  temperature  profile  and  has  less  effect  in  calcu¬ 
lating  the  wave  velocity.  However,  for  any  spot  size  which 
is  smaller  than  .5  cm,  the  conductive  heat  loss  through  radial 
direction  dominates  the  temperature  profile  and  thereby  affects 
the  wave  speed.  It  is  found  that  at  higher  absorbed  flux 

level  k„d>  >  1x10  watts/cm  and  large  beam  diameter  (R  >  0.5 
L  o 

cm)  a  simplified  linear  theory  or  phenomenological  analysis 
can  be  used  in  predicting  the  wave  velocity. 
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3.  A  NUMERICAL  STUDY  OF  SO.  IF,  POTENTIAL 

TWO-DI.  lEIJSIOHAL  AND  TARGET- INDUCED  IIYDRODYUAIIIC 
EFFECTS  Oil  LSC  WAVE  PROPAGATION 

3 .  1  IWTRODUCTIOM 

In  Section  2,  and  in  our  previous  studies,  wo  have  pre¬ 
sented  a  theoretical  formulation  which  seens  to  adequately 
reproduce  the  experimentally  observed  features  of  LSC  wave 
propagation.  In  particular,  it  is  shown  in  Figure  9,  tiiat 
our  model,  which  includes  radiative  pre-heating  as  the  pri¬ 
mary  transport  mechanism  of  the  wave  predicts  both  the  magni¬ 
tude  and  scaling  of  LSC  velocity  with  absorbed  flux  reasonably 
accurately.  Furthermore,  we  show  in  Figure  10,  that  the  vari¬ 
ation  of  wave  velocity  with  beam  radius  can  be  approximated 
as  nearly  linear  for  a  small  range  of  beam  radii  near  the  LSC 

row! 

throsnold  as  found  by  .  1SNW.  L  J  We  also  indicate  lower  theo¬ 
retical  predictions  for  the  LSC  velocity  than  experiment  at 
tiie  larger  beam  diameters.  Since  the  experimental  data  are 
sparse  and  the  beam  radius  variation  is  not  great,  the  transi¬ 
tion  to  independence  of  vq  on  beam  radius  lias  not  yet  been 
verified  experimentally.  As  has  been  stated,  this  model 
embodies  three  principal  simplifying  assumptions.  They  are 
(1)  the  model  is  one  dimensional,  (2)  the  deflagration 
approximation  is  used,  i.e.,  the  wave  is  in  pressure  balance 
with  tiie  atmosphere  at  p  =  1  atm,  and  (3)  the  wave  is  far 
from  the  target  so  that  wave  interactions  with  the  solid 
boundary  imposed  by  the  target  are  negligible. 

In  spite  of  the  previously  mentioned  reproduction  of 
the  experimental  observations  by  our  model,  t.iere  are  certain 
features  of  the  experimental  observations  which  are  not  pre- 

r  3  i 

dieted.  As  has  been  pointed  out  by  Klosterman  and  Byron,  J 
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a  s  Li.iplificd  one-dimensional  theory  produces  results  for  the 
wave  propagation  velocity  which  are  less  tuan  the  experimental 
t.ii Msureuents .  Including  all  of  the  non-linear  effects  and 
associated  loss  mechanisms ,  however,  we  find  the  theoretical 
predictions  of  velocity  to  be  in  reasonable  agreement  with 
tae  experinental  ouservations,  c.f..  Figure  9.  Since  a  1-D 
theory  should  place  an  upper  bound  on  the  LSC  velocity, 
Klosterman  and  Byron  argue  that  2-D  effects  must  be  respon¬ 
sible  for  their  observed  discrepancy.  The  small  discrepancy 

r  3 1 

in  itself  as  little  reason  for  concern,  but  as  has  been  shown, 
the  discrepancy  becomes  larger  as  the  laser  beam  diameter  is 
increased  -  in  apparent  disagreement  with  physical  intuition. 
Indeed,  Holographic  records  obtained  in  the  excell 'nt  study 
by  -15UU  indicate  that  the  wave  front  is  anything  but  planar, 
and  1-D,  and  that  the  LSC  wave  does  expand  radially  beyond 
the  diameter  irradiated  by  the  laser  beam.  Thus,  a  more 
detailed  investigation  into  the  2-D  effects  on  j_.SC  wave  propa¬ 
gation  is  in  order.  Also,  certain  other  features,  such  as 
target  effects  have  not  yet  been  studied. 

In  summary,  2-D  effects  can  be  produced  by  the  follow¬ 
ing  mechanisms: 

1.  A  focused  laser  beam  geometry,  wherein  the 
LSC  wave  propagates  from  regions  of  higher 
laser  fluence  into  regions  of  lower  laser 
f luence . 

2.  Due  to  non-perfect  attainment  of  pressure 
balance  with  the  atmosphere ,  on  the  time 
scale  of  the  interaction,  regions  of  higher 
than  atmospheric  pressure  can  exist  in  the 
heated  plasma.  Hence,  the  plasma  can  expand 
radially,  tiius  reducing  the  mass  of  gas 
heated  by  the  laser  beam. 

3.  The  incident  laser  pulse  is  not  necessarily 
spatially  uniform  as  a  function  of  radius. 

Typically  the  beam  is  Gaussian,  or  near- 
Gaussian.  Thus,  the  central  portions  of  the 
wave  are  subject  to  greater  heating  than  the 
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periphery,  resulting  in  radial  gradients  in 
the  thermodynamic  and  kinematic  variables. 

It  is  not  clear  to  us,  however ,  that  these  are  the  only 
effects  which  may  result  in  inadequacies  of  the  1-D  theories 
Lo  reproduce  all  details  of  the  experimental  observations. 

One  additional  feature  of  the  model  which  warrants  further 
investigation  is  tnat  associated  with  the  target-free  propa¬ 
gation  assumption.  Thus,  two  additional  complications  nay 
arise.  They  are: 

4.  As  the  ambient  gas  is  entrained  by  the  LSC 
wave,  it  is  accelerated  sharply.  Hear  a 
solid  surface,  the  heated  gas  is  accelerated 
into  the  surface,  stagnating  at  the  center- 
line  and  then  flowing  radially  out  across 
the  surface.  This  in  itself,  is  a  2-D  ef¬ 
fect.  However,  a  1-D  effect  results  as  well. 

It  is 

5.  A  compression  wave  which  propagates  back  onto 
the  target  surface  is  reflected  by  the  surface 
(since  zero  axial  velocity  must  obtain  at  the 
wall)  and  propagates  back  toward  the  propa¬ 
gating  CjSC  .  Hence,  an  LSC  which  lies  within 

a  round  trip  acoustic  transit  time  from  the 
target  is  subject  to  acceleration  by  the  re¬ 
flected  wave.  This  wave  can  propagate  at  a 
steady  velocity  and  thus  appears  independent 
of  the  target,  however,  the  velocity  of  such 
a  wave  can  be  as  much  as  an  order  of  magnitude 
higher  than  the  free-space  LSC  wave  subject  to 
the  same  laser  flux. 

While  a  detailed  study  of  each  of  these  effects  is  re¬ 
quired  before  firm  conclusions  regarding  their  relative  and 
absolute  importance  can  be  reached,  we  have  not  carried  out 
such  an  exhaustive  study  here.  A  full  2-D  treatment  is  re¬ 
quired.  Instead,  we  have  chosen  to  investigate  a  subset  of 
these  effects  from  which  pertinent  conclusions  can  be  reached 
by  use  of  a  quasi-  one-dinensional  code;  i.e.,  1-D,  but  with 
variable  laser  beam  area.  Hence,  we  address  in  the  following 
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sections ,  itens  (1),  (2)  and  (5).  Our  effort  is  directed  to 
an  understanding  of  the  propagation  of  LSC  waves  and  not  on 
tae  ignition  per  se .  Hence,  in  our  numerical  calculations 
as  discussed  in  Section  3.2,  we  assume  that  a  partially  ion¬ 
ized  plasma  which  has  boon  produced  by  energy  transport  from 
the  heated  target  exists  in  front  of  the  target.  The  results 
of  our  calculations,  which  are  presented  in  Section  3.2  and 
3.3  arc  based  on  modification  of  a  one-dimensional,  variable 
area  code,  HELVA,  developed  previously  by  Chapman  at  S3.^2(^ 
Our  modifications  to  this  code  for  the  present  study  are 
detailed  in  Section  3.4. 

3.2  TECHNICAL  DISCUSSION  Aih)  RESULTS 

Calculations  were  performed  on  a  laser  supported  com¬ 
bustion  wave  propagating  down  a  cylinder  or  conoidal  frustrum. 
TiU‘  shape  of  those  geometrical  volumes  is  determined  by  the 
characteristics  of  the  assumed  focusing  system  of  the  laser. 
The  initial  conditions  for  such  a  system  arc  shown  in  Figure 
11.  Here,  we  have  the  laser  shining  from  the  right  onto  the 
ambient  gas.  It  passes,  virtually  unattenuated,  through  the 
cold  zone,  and  is  absorbed  in  the  transition  and  hot  zones. 

The  hot  zone  is  initially  set  to  some  appropriately  high  tem¬ 
perature  (c.g.,  15,000  °K)  and  the  temperature  is  allowed  to 
decline  linearly  through  the  transition  zone  to  the  ambient 
temperature  of  300  °K.  The  length  of  the  transition  zone 
(typically,  one  centimeter)  and  the  initial  temperature  of 
the  hot  zone  was  not  found  to  be  very  critical  in  the  final 
results  obtained. 

One  of  the  critical  parameters  of  the  calculation  is 
the  initial  size  of  the  hot  plasma,  e.g. ,  the  ignition  mech¬ 
anism.  The  reason  for  this  sensitivity  is  due  to  the  compli¬ 
cated  gasdynamical  wave  interactions  which  will  now  be  dis¬ 
cussed  . 
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We  illustrate  the  evolution  and  propagation  from  an 
initially  hot  plasma  of  length  1  cm  in  Figure  12(a)-(d). 

Hero,  tiie  laser  is  incident  at  a  flux  of  3'105  watts/cri2 
onto  a  focal  spot  of  0.25  cm  radius.  The  laser  is  assumed 
to  lie  unfocused  here,  i.e.,  the  incident  rays  are  parallel. 

As  the  laser  is  turned  on,  the  plasma  begins  to  expand  by 
hydrodynamic  processes  and  the  JuSC  wave  front  begins  to 
form  by  radiation  transport  into  the  cold  gas  ahead.  The 
radiation  front  (or  inverse  bleaching  wave)  compresses  and 
accelerates  the  gas  ahead  of  it.  The  acceleration  leads  to 
the  formation  of  a  gasdynamic  shock.  Hence,  there  exists 
two  transition  regions  in  the  flow  field  -  the  gasdynamic 
precursor  shock  and  the  LSC  which  serves  as  a  "leaky  piston" 
and  supports  the  shock  wave.  Simultaneously  with  the  forma¬ 
tion  of  the  LSC  front  and  the  precursor  shoe];  wave ,  a  shock 
is  formed  at  the  front  surface  of  the  plasma  which  propagates 
back  toward  the  target  surface.  This  shod;  is  formed  as  a 
result  of  the  momentum  imparted  to  the  plasma  surface  by  the 
gas  which  expands  away  from  its  surface  due  to  the  deposition 
of  laser  radiation  and  subsequent  acceleration  away  from  the 
target  by  the  precursor  shock.  Since  the  flow  velocity  behind 
the  LSC  front  is  directed  towards  the  target,  the  shock  easily 
propagates  back  and  reflects  from  the  target  surface.  The 
shock  then  propagates  back  through  the  plasma,  but  cannot 
overtake  the  LSC  or  precursor!  shock  due  to  the  high  flow 

velocity  of  the  plasma  towards  the  target.  However,  the  exis- 

\ 

tonce  of  the  trailing  shock  servos  to  "tell  the  flow  ahead" 
of  the  presence  of  the  wall.  |  This  can  be  seen  in  Figure  12(d), 
where  the  LSC  velocity  increases  sharply  at  about  the  time 
the  secondary  shock  reflects  tarom  the  target  surface.  Two 
important  conclusions  are  appirent  from  Figure  12(d);  (1)  the 

reflected  shock  communicates  the  presence  of  a  target  (reflec¬ 
ting  surface)  to  the  LSC  wave  //hick  propagates  at  a  higher 
velocity  than  if  the  surface  wj  re  absent,  and  (2)  the  precur¬ 
sor  shock  propagates  at  constant  velocity  with  the  energy  re¬ 
quired  to  support  such  propagation  supplied  by  the  piston  (LSC) 


wave  system 
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vx  i  the  absorption  of  laser  radiation.  it  can  be  shown  that 
Lhe  Mechanical  energy  is  a  small  fraction  of  the  total  energy 
absorbed  by  the  plasma.  Furthermore,  the  presence  of  the  pre¬ 
cursor  Shock  produces  a  high  pressure,  high  density  region 
into  which  the  LSC  propagates.  The  formation  and  propagation 
of  the  LSC-shock  is  shown  in  Figure  12 (a) -(c)  where  we  plot 
temperature  and  pressure  as  a  function  of  distance  from  the 
target  for  times  of  22.02,  66.48  and  146.3  psec  respectively. 
Et  can  be  seen  that  the  LSC  wave  propagates  into  a  "modified 
ambient"  pressure  of  about  2.5  atm,  as  opposed  to  the  one 
atmosphere  assumption  implicit  in  the  common  deflagration 
assumption.  However,  it  can  be  seen  that  the  pressure  is 
maintained  nearly  constant  throughout  the  wave  which  is  con¬ 
sistent  wi til  the  deflagration  approximation. 


Before  leaving  Figure  12,  an  additional  feature  is  in 
eluded  on  Figure  12(d)  which  should  be  noted.  V/e  have  shown 
the  position  of  the  LSC  front  as  determined  arbitrarily  from 
the  location  of  the  8,000  «K  isotherm,  i.e.,  the  position  of 
maximum  brightness  to  the  visible  detector.  Also,  for  com¬ 
parison  we  show  the  position  of  a  trailing  compression  front 
waich  in  fact  corresponds  to  the  region  of  maximum  tempera¬ 
ture  behind  the  wave.  The  separation  between  the  two  corres¬ 
ponds  to  the  region  of  photon  absorption  and  thus  is  the 
absorption  length  of  the  laser  radiation. 


For  purposes  of  comparison,  we  include  in  Figure  13, 
an  analogous  calculation  at  a  laser  flux  of  5><105  watts/cn2. 
Due  to  the  larger  laser  flux,  it  can  be  seen  that  the  secon¬ 
dary  shock  reflects  from  the  target  in  a  shorter  time.  Also, 
we  include  on  Figure  13(c)  in  addition  to  the  temperature 
and  pressure,  the  density  profile.  The  existence  of  the  pre¬ 
cursor  shock  is  evidenced  by  the  sharp  rise  in  p  and  p  at 
about  9.6  cm.  The  pressure  remains  constant  at  about  3.5  atm 
through  the  LSC  wave  front,  whereas  the  density  decreases 
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sharply  as  the  temperature  increases  due  to  laser  heating. 
Evidence  of  subsequent  compression  of  the  plasma  is  apparent 
as  Lae  wave  cools  behind  the  maximum  temperature  region. 

A  summary  of  our  calculations  is  shown  in  figure  14, 
where  we  plot  LSC  wave  velocity  as  a  function  of  flux.  For 
comparison,  we  plot  the  experimental  correlation  curve  ob¬ 
tained  at  fISIJW;  v  =  13  Fpd  (cm/sec)  for  flux  density 

2 

h  (kw/cm  ),  pressure  p(atrn)  and  beam  diameter  ci  (cm)  for 
p  =  1  atra  and  d  =  0.5  cm. 

The  first  observation  is  that  our  computed  velocities 
are  considerably  higher  than  the  experimental  curve.  This  is 
due  directly  to  the  influence  of  the  target  on  our  calculated 
results.  As  we  have  indicated  in  Figures  12  and  13,  the  tar¬ 
get  influence  persists  to  distances  of  at  least  10  cm  and  a 
tine  scale  of  at  least  200  psec  after  LSC  formation.  The 
lGU'.J  measurements  were  made  at  similar  times  and  distances 
and  it  was  maintained  that  "target  free  propagation"  was 
observed  due  to  the  steady  propagation  of  the  LGC .  This  dis¬ 
crepancy  with  experiment  will  be  di„  ussed  in  Section  3.3. 

Secondly,  we  show  two  velocity  curves  entitled  long 
sone  and  short  zone.  The  long  zone  calculation  corresponds 
to  an  initial  plasma  length  of  7  cm,  whereas  the  second 
curve  corresponds  to  a  short  zone  of  1  cm  length.  The  long 
zone  curve  shows  a  lower  propagation  velocity  than  the  short 
zone  calculation.  This  result  indicates  the  sensitivity  of 
the  calculated  results  to  the  initial  (ignition)  conditions. 
The  long  zone  case  remains  more  free  from  the  target 
(although  not  entirely  free)  than  the  short  zone  due  to  the 
longer  time  for  the  secondary  shoe];  to  reflect  from  the  wall 
and  its  attenuation. 

Finally,  we  numerically  "force"  the  computer  code  to 
reproduce  the  deflagration  approximation  (of  constant  one 
atmosphere  propagation)  by  artificially  compelling  the  pres¬ 
sure  to  remaind  at  one  atmosphere.  This  technique  serves 
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liyure  14.  Plot  a  LSC  v/ave  velocity  vs  incident  laser  flux 
and  comparison  with  experimental  data. 
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Lwo  purposes;  (1)  by  maintaining  the  pressure  at  one  atmos¬ 
phere,  shoe);  formation  and  target  induced  hydrodynamic  pro¬ 
cesses  are  eliminated,  thus  the  target  free  propagation  mode 
is  reproduced,  and  (2)  mass  is  "artificially"  lost  from  the 
system,  albeit  at  a  rate  which  is  not  necessarily  real.  The 
results  obtained  by  this  method  of  solution  reproduce  the 
•ASiJU  correlation  rather  well.  Finally,  we  also  show  for  com¬ 
parison,  results  from  our  non-linear  solution  obtained  by 
the  quasilinearization  model  discussed  in  Section  2.  Both 
calculation  techniques  reproduce  the  results  identically. 

One  further  calculation  was  performed  in  which  we 
placed  a  constant  pressure  (1  a to)  boundary  condition  on  the 
left  end  of  the  "tube",  i.e.,  at  the  target  surface.  The 
secondary  wave  is  thus  not  allowed  to  reflect.  Results  ob¬ 
tained  by  this  technique  lie  midway  betv/een  the  "long  zone" 
and  "pressure  uniformly  equal  to  one  atmosphere"  results. 

In  order  to  convince  ourselves  that  the  initial  tem¬ 
perature  distribution  was  not  too  artificial  and  not  having 
an  exaggerated  effect  on  the  fluid  dynamics,  a  calculation 
was  performed  for  fifteen  microseconds  with  the  pressure 
lie  Id  to  one  atmosphere,  and  then  the  calculation  was  con¬ 
tinued  in  the  normal  way.  The  results  seem  to  be  quite 
similar  for  thermodynamic  profiles  and  for  the  LSC  wave  tra¬ 
jectory. 

3.3  DISCUGSIOII  OF  FOCUSED  GE0I1ETRIC  ADD  LATERAL  EXPANSION 

EFFECTS  — " 

As  was  indicated  in  Section  3.1,  two  potential  2-D 
effects  which  we  have  studied  are  expansion  due  to  focused 
laser  beam  geometry  and  radial  expansion  of  the  plasma 
beyond  the  confines  of  the  laser  channel. 
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If  the  beam  is  focused,  the  laser  flux  will  be  greater 
ihvu  Liie  focus  (say  at  the  target  surface)  than  at  some  dis¬ 
tance  away  from  the  target.  The  LSC  will  thus  propagate  into 
regions  of  lower  flux,  thus  explicitly  reducing  the  velocity. 
Fur thermore ,  due  to  the  nonparallel  channel  into  which  the 
LoC  propagates,  the  area  through  which  flow  occurs  increases. 
The  flow  process  corresponds  to  variable  area  flow  in  a 
channel  with  heat  addition.  Since  the  flow  is  subsonic,  the 
area  increase  reduces  the  flow  velocity.  Results  obtained 
froia  a  series  of  calculations  are  shown  in  Figure  14.  Here 
we  plotted  the  position  of  the  LSC  front  for  the  same  laser 
flux  at  the  focus  but  for  channels  with  divergence  half-angles 
of  0°,  2.58°,  and  5.43°  respectively.  As  can  be  seen,  the 
trajectories  are  very  similar  for  about  the  first  20  usee  of 
propagation.  After  that,  however,  the  LSC  velocity  becomes 
progressively  smaller  as  the  divergence  angle  is  increased. 

In  order  to  interpret  these  results,  we  must  separate 
the  effect  of  the  reduced  flux  at  a  particular  channel  loca¬ 
tion  from  the  purely  hydrodynamic  effect  of  increased  channel 
area.  If  v  «  Fr  where  F  is  the  laser  flux  and  r  the  channel 
radius,  we  find  that  simply  from  the  reduced  flux,  that  for 
a  divergence  half  angle  of  ^5°,  we  would  expect  a  velocity 
reduction  of  about  75%  from  the  velocity  of  constant  area 
flow.  Within  the  uncertainty  of  locating  the  position  of 
the  front,  we  find  from  Figure  15,  that  the  velocity  reduc¬ 
tion  is  about  the  same.  We  conclude,  therefore,  that  the 
udial  hydrodynamical  expansion  due  to  the  increased  channel 
area  has  negligible  influence  on  the  propagation  velocity. 

One  further  aspect  of  two-dimensional  flow  was  studied. 
That  is  the  radial  expansion  through  the  lateral  boundaries 
of  the  propagation  channel  defined  by  the  laser.  As  was  indi¬ 
cated  earlier  and  as  is  evident  from  Figures  12  and  13,  the 
pressure  behind  the  LSC  wave  and  between  the  precursor  shock 
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Figure  15.  Distance-time  plot  for  beam  divergence  angles  of 
0°,  2.577° ,  and  5.427°. 
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and  the  LSC  wave  is  considerably  greater  than  atmospheric, 
due  to  this  pressure  difference,  it  is  expected  that  flow 
v/ill  be  induced  in  the  radial  direction.  In  order  to  obtain 
a  reasonable  estimate  -  or  to  establish  an  upper  bound  on 
this  effect  -  we  chose  to  implement  the  model  which  is  pre¬ 
sented  in  Appendix  C.  We  reason  that  the  maximum  velocity 
at  which  the  flow  could  expand  rapidly  would  be  if  the  pres¬ 
sure  in  the  cylindrical  plasma  would  induce  a  cylindrical 
blast  wave  which  propagates  radially  away  from  the  plasma. 

The  model  presented  in  Appendix  C  was  used  to  calculate  the 
velocity  at  the  channel  boundary  induced  by  such  an  outwardly 
propagating  blast.  This  velocity  is  then  used  to  compute  the 
mass  convected  through  the  lateral  channel  boundary,  which  is 
thus  removed  from  each  computational  zone  during  the  calcula¬ 
tion.  This  model  was  used  for  all  of  our  calculations,  thus 
its  effect  is  included  in  all  of  the  results  included  herein, 
o.cj..  Figures  12,  13,  14.  Accordingly,  we  find  that  even 
including  such  radial  mass  loss  from  the  flow  field,  our 
calculations  do  not  agree  with  the  MSNW  correlation  (except 
for  the  p  =  1  atm  calculation) .  This  indicates  that  the 
radial  mass  loss  is  not  adequate  to  reduce  the  pressure 
behind  the  precursor  shock-LSC  system  over  the  distance  and 
time  scale  for  which  the  calculations  were  performed;  i.e., 

,,J  200  Usee  and  10  cm  from  the  target.  That  this  is  true  can 
be  seen  from  Figures  12  and  13  which  indicate  that  the  pres¬ 
sure  levels  remain  high  over  the  duration  of  the  calculation. 
For  instance,  we  find  that  for  the  case  shown  in  Figure  3 
(5*10J  watts/ern^,  R  =  .25  cm),  the  mass  loss  rate  from  a 
zone  between  the  shock  and  the  LSC  is  of  order  0.5  grams/sec. 
Since  the  mass  per  zone  is  about  10x10  grams,  the  charac¬ 
teristic  decay  time,  10*10-5  grams/0.5  grams/sec  =  2*10-4 
sec,  is  about  200  Msec.  That  is,  about  200  usee  subsequent 
to  passage  of  the  wave.  We  find  no  evidence  of  a  signifi¬ 
cant  pressure  reduction  behind  the  wave  in  our  calculations, 
however.  But,  we  do  find  that  the  combined  effect  of  a 
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i  ili.il  velocity  and  radial  conduction  is  of  suflioient:  mayni- 
Ludo  to  produce  a  plasma  extending  beyond  the  confines  of  the 
light  channel  as  found  in  the  holograms  and  framing  photo¬ 
graphs  obtained  by  MSNW. 

We  are  thus  faced  with  the  conclusion  that  according 
to  our  numerical  calculations,  the  2-D  effects  of  radial 
flow  and  channel  expansion  due  to  the  focal  geometry  are  not 
adequate  to  explain  the  experimental  results  obtained  by 
.15dW.  However ,  use  of  the  deflagration  approximation  brings 
our  calculations  into  good  agreement,  thus  indicating  that 
any  target  effects  present  in  the  ilSlh/  experiments  are  small 
indeed.  We  must  therefore  conclude  that  according  to  the 
experimental  data,  more  target  decoupling  occurs  in  reality 
than  is  predicted  by  a  one-dimensional  calculation.  The  two- 
dimensional  expansion  over  the  target  surface  must  weaken  the 
reflected  shock  wave  and  result  in  additional  decoupling  over 
the  time  scale  of  the  experiment. 

In  order  to  model  this  phenomena,  a  2— D  calculation  is 
desirable.  However,  we  feel  that  the  stagnation  point  flow 
field  is  undoubtedly  highly  complicated  by  turbulent  mixing 
with  the  expanding  target  material.  A  calculation  of  this 
complicated  phenomena,  including  a  reasonable  description  for 
the  target  vapor-shock  wave  interaction  is  beyond  the  scope 
oj  the  present  study.  Further  work  in  this  area  is  necessarv 
lief  ore  such  a  computation  will  bear  any  resemblance  to  reality 
and  thus  be  worth  pursuing. 

At  the  present  time,  however,  for  the  purpose  of 
systems  studies, interpretation  of  data  and  experimental  guid¬ 
ance,  we  feel  that  either  our  non-linear  model,  or  the  one- 
dinensional  variable  area  model  described  herein  when  coupled 
with  the  deflagration  approximation  is  adequate.  Inclusion  of 
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non-linear  effects  and  loss  terns  renove  nost  of  the  discre¬ 
pancies  noted  by  MSIJW.  The  results  obtained  by  us  as  dis¬ 
cussed  in  this  section  and  in  Section  2  seen  to  be  physically 
consistent,  and  when  interpreted  properly  are  consistent  with 
experimental  observations. 

One  area  which  has  not  been  studied,  however,  and  which 

could  be  addressed  by  technique  such  as  those  described  in 

Section  3,  is  the  transition  from  LSC  to  the  LSD  node.  This 

could  be  very  important  for  systems  implications  when  using 

repetitively  pulsed  lasers  operating  at  low  cw,  but  high- 

peak  powers.  We  have  shown  the  existence  of  a  precursor 

shock  whose  strength  increases  with  increasing  laser  fluence. 

The  electron-cascade  breakdown  behind  the  precursor  shock  is 

responsible  for  absorption  in  the  shock  zone  and  ultimate 

decoupling  of  the  shock  from  the  LSC.  The  flux  region  above 
G  2 

about  5-^10  watts/cn  will  thus  be  dominated  by  this  transi¬ 
tion  regime  which  evidences  combined  LSC  and  LSD  phenomena, 
and  must  be  studied  by  numerical  techniques  due  to  the  com¬ 
plexity  of  the  hydrodynamic  processes  involved. 


3.4  THE  COMBUSTION  WAVE  I10DEL 

The  problem  was  solved  from  the  numerical  integration 
of  the  Eulerian  conservation  equations  for  mass, 


3p 

Tt 


+  V* (pu)  =  0 


(28) 


for  momentum, 
J  ( pu) 


dt 


+  V(pu^)  +  Vp  =  0 


(29) 


and  for  energy. 


Jyt~G')'  +  +  P^’u  +  V*q 


=  0 


(30) 
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where  p  is  the  density,  u  is  the  velocity,  p  is  the  pressure, 
e  is  the  specific  internal  energy  and  cj  is  the  generalized 
luvi t  flux  (from  radiative  absorption  and  conductive  and  radia¬ 
tive  conduction) . 

To  integrate  these  terms  numerically,  an  integration 
scheme  developed  by  Chapman ^20]  has  been  used.  The  scheme  used 
is  explicitly  given  by  Chapman  and  differs  in  the  integration 
of  his  Equation  3.21: 


/;in+l  n+1  .n  n  \ 

Lii.Vj+’i  ~  -jh  V*) .  c  + »  +  i 


(31) 


whore  C  is  the  conduction  term, 

C  =  (PU  A)f^  -  (pu  A)n_!i 
e  j  e  j 

W  is  the  work  tern, 


(32) 


W  = 


-  (*j+i <i  -  vr )  ■ 

and  L  is  the  term  peculiar  to  the  LSC  problem, 


1  JR  +  JC  +  CR  +  CC  +  A 


(33) 


(34) 


Here  JR  and  are  the  leakage  of  energy  from  the  section  by 
lateral  radiative  and  conductive  diffusion,  CR  and  C are  the 
transmission  of  energy  between  mesh  cells  by  radiative  and 
conductive  diffusion, 


CR  = 


|K>;:ivi  +  <vr  *,][*& -*r]  (3s 

vJbr1  -  T"-i]  i/<2«>  - 


r  i  \n+l  .  ,  , 

l'Xr  j  Ai  +  ( 
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CC  = 


>  n+1 

C;j+1  Aj+1 

+  ‘V?"  A^ 

-  *r] 

>r  Ai  - 

vJI 

ipn+i  _ 

-  j 

(36) 


and  A  is  the  absorption  of  energy  from  the  laser  beam 


A  =  1-exp  -k"+1  Ax  AFn^  a. A1 
3  1+1  3+1 


(37) 


Fj+i  is  the  laser  flux  impinging  on  mesh  cell  j  from  the 
direction  of  j+1.  The  flux  is  attenuated  by  the  formula 


^n+1 

j 


Ax)]  Aj+i/Aj 


(38) 


whore  Fj+^  F0»  the  incident  flux.  The  notation  used  in 
defining  C  and  W  are  given  by  Chapman  and  the  formulas  for 


the  terms  contained  in  L  are  given  in  Section  3.2  of  Boni, 
Cohen  and  Su. 


The  energy  difference  equation  is  solved  iteratively 

since  W  and  L  are  functions  of  e1?*,1.  This  is  due  to  the 

J  *  *5 

fact  that  the  work  term  contains  pressure, 


pn+i  _  p/nn+l  n+1  \ 

pj+>i  p(Vj  '  ej+d  - 

and  L  is  a  functior  of  the  temperature, 


Tn+1  =  T  /nn+i  Pn+1  \ 
3  +  ’i  T(Pj+ta  '  ej+li) 


n+1 


n+1 


Those  two  intrinsic  functions  are  determined,  in  the  case  of 
air,  by  a  simple  fit  to  Predvoditelev ' s  data,^1^ 

e  =  5.036^±03  T2  -  P/p  , 
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4.  EXPERIMENTAL  INVESTIGATION  OF  THE  EFFECTS  OF  PART’Tr’nr  nmn 

matter  on  short  wavelength  pulsed  laser  propagation 

4 ■ 1  INTRODUCTION  and  background 

[,21  Durin9  the  past  several  years,  theoretical  work  at 

s  as  WGl1  as  work  at  the  Air  Force  Weapons  Laboratory123’ 

has  been  reasonably  successful  as  a  framework  for  the  inter¬ 
pretation  of  observations  based  on  dirty-air  breakdown  experi¬ 
ments  conducted  at  Avco-Everett  Research  Laboratory'24’  Lin¬ 
coln  Laboratory'  and  at  United  Aircraft  Research  Labora¬ 
tories  .  Theoretical  calculations  have  been  performed 
which  are  based  on  a  model  of  particle  induced  breakdown  which 
considers  vaporization  of  the  particle  due  to  laser  heating, 
followed  by  an  electron  cascade  process  initiated  either  in 
the  heated  particle  vapor  or  in  the  adjacent  air  which  has 
been  shock-heated  as  a  result  of  the  rapidly  expanding  particle 
vapor .  It  has  been  pointed  out  that  this  picture  of  the  dirty- 
air  breakdown  problem  is  similar  in  many  respects  to  the  target- 
interaction  problem  except  for  the  different  geometry  and 
spitial  scale.  Indeed,  the  experimental  LSD  thresholds  as 
observed  by  Hall  for  plane-target  interactions  -  "infinite 
diameter  particles"  -  are  about  the  same  as  those  observed  in 
Lac  above  mentioned  air  breakdown  experiments.  Specifically 
a:i  we  have  pointed  out  previously'26’  the  laser  fluenco  associ¬ 
ated  with  containment  breakdown  in  air  and  for  the  LSD  thresh¬ 
olds  are  both  of  order  10  joules/cm2  for  pulse  times  of  order 
10  sec  and  are  consistent  with  an  infinite  diameter  particle 
extrapolation.  In  any  event,  it  has  been  shown  both  theoreti¬ 
cally  and  experimentally  that  the  breakdown  threshold  for  a 
particle-laden  atmosphere  may  be  several  orders  of  magnitude 
less  than  the  clean  air  breakdown  threshold  of  about  3  x  10s 
watts/cm2  for  C02  wavelengths  depending  on  pulse  length  and 
the  particulate  concentration  and  particle  size  distribution. 
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The  main  thrust  of  the  effort  in  the  community  over 
the  past  several  years  has  been  to  correlate  and  to  under¬ 
stand  both  the  long-pulse  and  short-pulse  experimental  data 
for  CO 2  laser  propagation,  i.e.,  at  10.6|j  including  the  effect 
of  atmospheric  particulate  impurities.  However,  since  it  has 
been  shown  that  high  fluence  C02  laser  propagation  through 
the  atmosphere  is  limited  due  to  non-linear  effects,  viz, 
thermal  blooming  and  gas  breakdown,  and  that  more  optimum 
propagation  characteristics  may  exist  at  reduced  wavelengths, 
a  large  amount  of  work  has  been  initiated  to  develop  a  high- 
intensity  laser  source  which  operates  efficiently  at  wave¬ 
lengths  lower  than  10. 6y.  Accordingly,  a  significant  amount 
of  work  has  gone  into  developing  chemical  lasers,  such  as  the 
DF  laser,  cold  cathode  electron  beam  stimulated  lasers  such 
as  the  CO  laser,  hybrid  chemical-electrical  excitation  lasers, 
and  finally  most  recently,  relativistic  electron  beam  stimu¬ 
lated  lasers  which  operate  in  the  visible  and  UV  regions  of 
the  spectrum. 

Based  on  inverse  bremstrahlung  scaling  as  discussed  in 
Appendix  D,  it  is  anticipated  that  breakdown  thresholds  at 
reduced  wavelength  should  scale  directly  with  the  square  of 
the  laser  frequency  w.  Accordingly,  the  non-linear  problem 
associated  with  gas  breakdown  may  be  less  restrictive  at  lower 
wavelengths  and  a  higher  fluence  laser  pulse  :aay  be  propagated. 
Furthermore,  it  is  expected  that  the  breakdown  threshold  flux 
increases  as  the  laser  pulse  length  becomes  shorter.  These 
observations  follow  from  the  classical  limit  of  a  simple 
energy  balance  equation,  c.f.  Appendix  D,  which  relates  the 
rate  of  growth  of  average  electron  energy  to  the  electromag¬ 
netic  field  of  the  laser. 

It  is  obvious  that  many  potentially  complicated  effects 
could  and  have  been  included  in  calculations  such  as  those 
carried  out  by  us  previously  for  CO2  laser  interactions.  How¬ 
ever,  the  utility  of  a  simplified  approach  is  highly  desirable. 
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The  importance  of  timely,  accurate  and  relevant  experimental 
data  in  the  guidance  and  formulation  of  theoretical  model 
development  is  recognized  as  essential  for  the  attainment  of 
a  reliable  solution  to  the  objectives  of  the  DARPA  laser  pro¬ 
gram.  For  this  reason,  and  because  of  the  importance  of  ob¬ 
taining  answers  to  the  wavelength  and  pulse  length  scaling 
question,  we  have  performed  a  small  number  of  measurements 
related  to  dirty-air  breakdown  at  Nd- laser  wavelengths .  The 
choice  of  the  Nd  laser  ( 1 . 0 6 y  wavelength)  was  based  on  the 
availability  of  an  operational,  high  fluence  laser  facility 
at  reduced  wavelength  and  pulse  length  so  that  interim  informa¬ 
tion  could  be  obtained  prior  to  the  availability  of  appropriate 
gas  laser  systems. 

Experimental  data  are  obtained  and  results  presented 
for  breakdown  of  filtered  air,  Si02  fibers  ( 9 u )  and  carbon 
fibers  (6.6p  and  68p)  by  30  nsec,  1.06y  laser  pulses.  It  is 
shown  that  the  results  correlate  well  with  data  obtained  at 
10.6m  and  that  the  breakdown  thresholds  scale,  within  the 
experimental  uncertainties,  with  the  square  of  the  laser  fre¬ 
quency  between  1.06p  and  10. 6p.  The  lowest  threshold  measured 
is  for  the  68|i  carbon  fiber  and  is  2.3  x  10 10  watts/cm2  for 
the  30  nsec  pulse  length.  Lower  values  are  expected  for 
longer  pulses.  The  smaller  diameter  impurities  result  in  in¬ 
creased  breakdown  thresholds  approaching  the  flux  limited  value 
of  3  x  10 1 1  watts/cm2  for  clean  air.  As  has  been  found  for 
breakdown  at  10. 6y ,  we  find  an  independence  of  breakdown 
threshold  on  the  material  composition  of  the  impurity. 
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4 . 2  EXPERIMENTAL 

The  test  matrix  included  a  limited  number  of  experiments 
related  to  dirty-air  environments  but  was  principally  concerned 
with  breakdown  of  single,  small  diameter  fibers.  We  performed 
measurements  on  Si02  fibers  (9y  diameter)  and  carbon  fibers 
(6.6|i  and  68p  diameter)  in  both  air  and  argon  atmospheres.  In 
addition,  the  breakdown  threshold  for  filtered  air  (O.ly  fil¬ 
ters)  was  measured.  As  mentioned,  the  experiments  were  con¬ 
ducted  by  S3  at  the  Battelle  Columbus  Laboratories  on  the  CGE 
640  Ud:  Glass  facility  with  the  laser  operated  in  the  single 
pulse,  30  nsec  FWHM  mode.  This  facility  is  capable  of  pro¬ 
ducing  up  to  450  joules  at  30  nsec,  although  we  used  signifi¬ 
cantly  less  energy  in  our  experiments.  A  schematic  layout  of 
the  experimental  arrangement  is  shown  in  Figure  17. 

A  one  foot  diameter,  one  foot  high,  evacuable  metal 
cylinder  was  used  as  the  test  chamber.  The  chamber  was 
equipped  with  a  number  of  ports  around  its  circumference  for 
diagnostic  access.  The  main  laser  pulse  is  brought  into  the 
chamber  unfocused,  with  a  diameter  of  7  cms .  All  of  the  beam 
focusing  was  performed  inside  the  chamber  with  a  lens  of  focal 
length  3  meters.  An  extension  port  to  facilitate  the  long 
focal  length  lens  was  provided.  The  flux  density  at  the  focal 
volume  was  varied  by  changing  the  laser  operating  parameters, 
i.e.,  by  varying  the  energy  input  and  number  of  amplification 
stages,  rather  than  by  repeated  changes  of  the  focusing  lens 
and/or  the  size  and  location  of  the  focal  volumes  which  v/as 
maintained  constant  throughout  the  experiment.  Laser  beam 
diagnostics  include  incident  energy,  beam  size  and  spatial 
distribution,  and  time  dependence  of  the  incident  pulse.  The 
incident  energy  and  pulse  time  dependence  and  focal  position 
v/as  monitored  for  each  shot. 

Diagnostics  of  the  breakdown  event  for  each  shot  in¬ 
cluded  the  following:  (1)  a  measurement  of  the  transmission 
of  the  pulse  through  the  focal  region,  and  the  time  dependence 
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of  the  transmitted  pulse,  (2)  a  measurement  of  che  time- 
resolved  emission  from  the  breakdown  region,  (3)  open-shutter 
photography  of  the  breakdown  region,  (4)  streak  camera  photo¬ 
graphy  of  the  breakdown  region,  and  (5)  time  resolved  cross 
beam  absorption  by  use  of  a  He-Ne  laser. 

It  was  determined  that  the  transmitted  laser  pulse 
and  plasma  emission  measurements  were  good  indicators  of  the 
breakdown  event.  The  streak  camera  records  were  useful  in 
establishing  the  location  of  plasma  luminosity  at  distances 
away  from  the  fibers;  however,  a  time-delay  fault  in  the 
camera  (of  123  nsec)  did  not  allow  the  plasma  motion  during 
t.he  laser  pulse  to  be  measured.  Since  the  program  funding 
provided  for  only  a  one  week  stay  at  Battelle,  the  problem 
was  not  resolved  and  accordingly  the  streak  camera  records 
were  of  limited  usefulness.  One  additional  difficulty  encoun- 
'orod  with  the  experiment  was  with  the  cross-beam  transmission 
measurement.  It  was  discovered  that  due  to  a  predominance  of 
fiber  material  in  the  plasma,  plasma  emission  at  6328A  was 
sufficiently  great  that  meaningful  absorption  measurements 
could  not  be  obtained.  In  retrospect,  a  polarizer  placed  in 
front  of  the  photodiode  should  be  adequate  to  eliminate  the 
unpolarized  plasma  emission  from  interference  with  the  polar¬ 
ized  He-Ne  radiation. 

In  any  event,  good  and  reproducible  data  were  obtained 
on  laser  transmission  through  the  focal  volume,  plasma  emis¬ 
sion  and  open  shutter  photography.  By  a  combination  and  inter¬ 
comparison  of  these  diagnostics  we  could  determine  whether 
breakdown  occurred  and  at  what  time  during  the  pulse,  the  per 
cent  transmission  through  the  focal  region  and  the  shape  of 
the  breakdown  plasma  so  produced.  Also,  measurements  of  the 
plasma  emission  have  been  found  to  be  useful  indicators  of 
breakdown  and  correlate  well  with  the  through-transmission 
measurements.  Due  to  the  emission  from  the  heated  and  vapor¬ 
izing  fiber,  the  emission  level  increases  during  the  laser 
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pulse  and  then  decays  well  after  the  laser  is  turned  off. 
However,  the  initiation  of  gas  breakdown  causes  a  significant 
increase  in  the  emission  level  over  and  above  that  present 
due  to  fiber  heating  and  vaporization.  The  duration  of  the 
"plasma  flash"  then  indicates  both  the  existence  and  the 
duration  of  the  breakdown  plasma. 

The  size  of  the  laser  beam  at  the  focal  point  was 
measured  by  calibrated  "burn  paper."  It  was  determined  that 
the  focal  spot  was  0.075  cm  in  diameter  (mean  diameter  for 
a  slightly  oblong  spot  of  .070  x  .080  cm)  at  the  one-half 
power  point. 

4 . 3  EXPERIMENTAL  RESULTS 

First  of  all,  we  discuss  features  of  our  observations. 
For  the  air  breakdown  experiments,  a  typical  plot  of  the 
incident  and  transmitted  laser  pulse  and  the  plasma  luminosity 
trace  is  shown  in  Figure  18.  It  is  shown  that  breakdown 
occurs  very  near  the  peak  power  point  as  evidenced  by  the  sharp 
reduction  in  the  transmitted  pulse  trace.  The  plasma  luminos¬ 
ity  trace  indicates  that  the  emission  increases  to  a  maximum 
and  then  decays  over  a  time  period  of  about  1  psec.  This 
value  of  decay  time  is  typical  of  all  our  results.  The  photo¬ 
diode  used  to  monitor  the  plasma  emission  was  displayed  on  a 
much  longer  time  scale  than  the  laser  pulse,  thus  the  emission 
rise  time  was  not  determined.  The  maximum  value  of  the  plasma 
intensity  and  the  decay  time  were  the  primary  parameters  of 
interest. 

Two  separate  methods  were  utilized  to  establish  the 
breakdown  threshold  and  will  be  illustrated  for  the  case  of 
the  filtered  air.  First  of  all  we  plotted  the  percent  trans¬ 
mission  thre  ’gh  the  focal  region  as  a  function  of  energy 
incident  onto  the  focal  spot.  As  is  shown  in  Figure  19,  the 
transmission  remains  at  a  constant  high  value  until  the 
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IS  joule  point,  at  which  a  sharp  break  occurs  and  the  trans¬ 
mission  drops  rapidly  to  less  than  10?,.  Also  on  Figure  19,  we 
plot  the  emission  from  the  breakdown  region  as  measured  by  a 
monitor  photodiode.  It  is  seen  that  below  15  joules,  no  de¬ 
tectable  radiation  is  emitted,  whereas  for  laser  energies 
greater  than  15  joules,  the  emission  from  the  plasma  rises 
rapidly  as  the  laser  fluence  is  increased.  This  is  the  visible 
spark  or  flash  which  occurs  upon  breakdown.  It  is  believed 
that  these  two  methods  utilized  simultaneously  (along  with  the 
observation  of  a  sudden  drop  in  the  transmitted  laser  pulse 
trace)  present  a  clear  and  unambiguous  determination  of  the 
existence  of  breakdown. 

The  energy  of  15  joules  onto  a  focal  spot  area  of 

4.4  x  10-3  cm2  for  a  30  nsec  pulse  yields  a  breakdown  threshold 

of  1.1  x  1011  watts/cm2.  Since  the  focal  spot  diameter  is 

about  0.075  cm,  it  is  believed  on  the  basis  of  theoretical  cal- 
[22] 

culations  that  diffusion  losses  are  negligible.  The  value 

for  breakdown  threshold  is  about  a  factor  of  2-3  lower  than 
the  commonly  accepted  value  for  clean  air.  Even  though  we 
filtered  the  inlet  air  by  use  of  a  O.ly  filter  and  evacuated 
the  test  chamber  before  admitting  the  test  gas,  the  reduction 
in  threshold  from  the  clean  air  value  indicates  that  presence 
of  some  impurity  particles  of  diameter  as  great  as  ly .  This 
observation  is  not  surprising,  if  one  recognizes  that  the  focal 
volume  subject  to  the  breakdown  fluence  in  our  experiments  was 
relatively  large  due  to  the  use  of  a  three  meter  focal  length 
lens.  Thus  the  peak  intensity  region  was  over  1  cm  long  and 
the  presence  of  a  very  low  number  density  of  impurities  could 
lead  to  the  existence  of  containment  particles  within  such  a 
large  focal  volume.  For  instance,  with  a  focal  volume  of 
4  x  10” 3  cm3,  a  particle  number  density  of  about  2.2  x  102  cm- 3 
would  lead  to  the  existance  of  one  particle  within  the  focal 
volume.  Indeed,  the  open  shutter  photographs  indicated  that  the 
breakdown  region  was  not  localized  to  one  spot  from  shot  to 
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shot  for  the  air  shots,  indicating  the  presence  of  particles 
within  the  focal  volume.  Also  we  saw  cases  of  more  than  one 
breakdown  within  the  focal  volume.  (Placement  of  the  fiber 
w 1  thin  the  focal  volume,  however,  always  localized  the  break¬ 
down  origin  at  the  fiber.)  In  any  event,  the  purpose  of  mea¬ 
suring  the  air  breakdown  threshold  in  these  experiments  was 
to  characterize  the  environment  in  which  the  fiber  breakdown 
experiments  would  take  place.  It  has  been  shown  that  the  back¬ 
ground  gas  breakdown  threshold  is  considerably  greater  than 
that  measured  for  the  smallest  fiber  and  thus  does  not  affect 
determination  of  the  particle  induced  breakdown  thresholds. 

Subsequent  to  characterizing  the  background  gas  environ¬ 
ment,  we  proceeded  to  determine  the  breakdown  threshold  of  the 
oi02  and  carbon  fibers.  A  data  reduction  technique  similar  to 
that  described  previously  for  the  filtered  air  was  employed. 

.  Iio  only  significant  difference  between  the  data  reduction  pro¬ 
cedure  for  the  fibers  and  that  for  the  air  is  that  the  fibers 
affected  the  nature  of  the  plasma  emission  curves  duo  to  the 
heating  and  vaporization  of  the  surface  of  the  fibers.  In¬ 
stead  of  a  non-de tec table  level  of  plasma  emission  until  break¬ 
down  occurred,  it  was  noted  that  the  plasma  emission  rose  slow¬ 
ly  with  increasing  laser  fluence.  Then,  at  the  threshold,  a 
break  in  the  curve  occurred  and  the  emission  rose  much  more 
rapidly  with  increasing  laser  fluence. 

One  otaer  interesting  feature  of  the  transmission  curves 
is  illustrated  in  Figure  20  for  the  case  of  the  63m  carbon 
iber.  Here,  we  present  two  transmitted  pulse  curves,  illus¬ 
trating  the  behavior  very  near  the  breakdown  threshold.  It 
is  shown  that  just  under  threshold,  breakdown  begins  at  the 
!>e  d:  power  point.  However,  the  breakdown  is  not  self-sustain¬ 
ing  and  dissipates  quickly  as  the  laser  fluence  decreases. 

Vhe  plasma  emission  remains  low  as  is  shown. 

Just  slightly  over  threshold,  however,  the  breakdown 
begins  before  the  peak  power  point  is  attained.  High  beam 
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attenuation  is  maintained  throughout  the  laser  pulse.  In  the 
example  shown,  the  threshold  is  not  exceeded  by  much,  however, 
and  partial  transmission  is  recovered  as  the  laser  fluence 
decreases.  Also  note  the  significant  increase  in  plasma  emis¬ 
sion  in  the  "just  over  threshold"  case. 

A  further  increase  in  laser  fluence  produces  full  break¬ 
down  and  the  transmitted  pulse  remains  near  zero  for  the  full 
duration  of  the  laser  pulse.  Near  zero  transmission  implies 
full  blocking  of  the  beam,  which  indicates  that  the  breakdown 
plasma  expands  rapidly  to  fill  the  cross  section  of  the  focal 
volume  and  subsequently  propagates  up  the  beam  as  a  LSD  wave. 

The  open  shutter  photographs  show  a  bright  luminous  column  of 
diameter  equal  to  the  focal  diameter  propagating  up  the  beam 
towards  the  laser  source. 

In  any  event,  we  have  established  the  breakdown  thresh¬ 
old  for  /58]i  carbon,  6.6p  carbon  and  9y  Si02  fibers  as  2.3  x  1010 
watts/cm2,  4.6  x  10 1 0  watts/cm2  and  6.1  x  10 6  watts/cm2  respec¬ 
tively  for  the  fibers  placed  in  air. 

The  carbon  fiber  experiments  were  repeated  in  a  fil¬ 
tered  argon  environment  also.  In  plotting  the  %  transmission 
and  plasma  emission  curves  versus  incident  laser  fluence,  it 
was  not  possible  to  distinguish  any  difference  for  the  break¬ 
down  thresholds  between  the  fibers  located  in  air  or  in  argon. 
~his  result  is  consistent  with  previous  measurements^28^  where 
it  was  found  t^at  the  breakdown  threshold  in  N2  and  Ar  were 
practically  identical  for  a  background  pressure  of  one  atmo¬ 
sphere  . 

In  Figure  21,  we  summarize  our  results  and  compare  with 
t/iose  of  other  investigators.  First  of  all,  we  use  the  results 
of  Lencioni  who  correlated  the  available  data  at  10. 6y  to 
obtain  a  plot  of  breakdown  flux  versus  laser  pulse  length,  in¬ 
cluding  the  dependence  on  particle  size.  In  order  to  compare 
with  the  1.06p  data,  we  have  scaled  the  10. 6y  curves  by  a 
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factor  of  100.  The  curve  shows  two  distinct  regimes  corres¬ 
ponding  to  short  and  long  pulse  propagation.  For  short  pulse 
propagation,  i.e.,  £  10“ 9  sec,  the  constant  energy  regime 
exists.  The  curve  shown  has  a  slope  of  1000  joules/cm2.  For 
pulse  times  greater  than  1  nsec,  the  curve  breaks  and  makes  a 
transition  to  the  constant  flux  regime.  It  is  shown  that  the 
effect  of  increasing  particle  diameter  is  to  decrease  the 
threshold  and  to  increase  the  pulse  duration  at  which  transi¬ 
tion  to  the  constant  flux  regime  occurs. 

The  data  obtained  by  us  is  plotted  at  the  30  nsec  pulse 
duration.  As  can  be  seen,  within  the  experimental  error,  which 
is  estimated  to  be  about  ±30%,  our  results  indicate  that  the 
1/A 2  scaling  according  to  the  inverse  bremsstrahlung  arguments 
seems  to  be  valid.  In  addition,  results  recently  obtained  by 
Lencioni^^  (after  our  experiments  were  completed)  indicate 
that  1/A2  scaling  is  valid.  One  very  important  feature  to 
note  here  is  that  in  addition  to  the  fact  that  both  our  work 
and  the  work  of  Lencioni  verifies  the  1/A2  scaling,  pulse 

length  scaling  information  can  be  obtained  as  well.  Once 
again,  then,  we  have  indication  that  the  breakdown  threshold 
versus  pulse  length  curves  are  similar  for  both  l.OCy  and  10. 6y 
with  the  appropriate  wavelength  scaling  included. 

Also  in  Figure  21,  we  have  included  for  completeness 

(and  to  give  some  credibility  for  extending  the  constant  energy 

curve  to  10  picosec)  the  results  of  two  recent  works.  First 

f  31] 

of  all,  the  work  of  Ireland,  et  al.  at  20  psec  is  in  ex¬ 

cellent  agreement  with  the  1000  j/cm2  extrapolation.  The  work 

r  32  ] 

of  McPherson  and  Gravel  at  3  nsec  is  slightly  lower  than 

the  1000  j/cm2  curve.  This  small  discrepancy  is  easily  ex¬ 
plainable,  however.  McPherson  and  Gravel  conducted  their  ex¬ 
periments  with  a  mode  locked  CO2  TEA  laser.  The  breakdown  did 
not  occur  until  the  5th  or  6+_h  pulse  in  the  train.  ""  'y  main¬ 
tain  that  there  is  evidence  that  the  first  pulses  of  the  train, 
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winle  not  causing  visible  breakdown,  may  have  caused  some  pre¬ 
ion  ization  of  the  air  in  the  focal  region  thus  lowering  the 
threshold . 

Finally,  in  Figure  22,  the  results  presented  in  Figure 
21  are  plotted  as  breakdown  threshold  as  a  function  of  particle 
diameter.  Again,  we  scale  the  10. 6y  results  by  a  factor  of  100 
(corresponding  to  the  30  nsec  pulse  width)  which  are  shown  as 
the  dashed  curve. 

In  summary,  we  can  draw  the  following  conclusions  from 
our  work: 

(1)  There  is  a  particle  size  effect  on  lowering  the 
breakdown  threshold  for  particles  between  ly 
(estimated  size  of  the  impurities  present  in  the 
filtered  air  experiments)  and  68y.  The  smaller 
particles  require  a  larger  investment  of  energy 
since  the  rapid  hydrodynamic  expansion  for  the 
smaller  particles  imposes  greater  energy  loss 
mechanisms.  The  larger  particles  expand  more 
slowly  and  represent  a  smaller  loss. 

(2)  At  least  for  the  two  materials  which  we  have  mea¬ 
sured,  there  is  an  independence  on  material  com¬ 
position  of  the  fibers. 

(3)  For  1  atm  ambient  pressure,  we  could  furthermore 
make  no  distinction  between  air  and  argon  as  the 
background  gas. 

(4)  The  breakdown  phenomena  are  similar  at  both  1.06y 
and  10. 6y  and  scale  as  1/A2  indicating  that  in¬ 
verse  bremsstrahlung  is  the  dominant  absorption 
mechanism. 

(5)  Comparison  of  our  results  with  those  of  Lencioni ^  . 

McPherson  and  Gravel [32]  and  with  Ireland,  et  al.l31J, 
indicate  that  pulse  length  scaling  is  similar  at 
both  1.06y  and  10. 6y. 

(6)  According  to  the  scaling  arguments,  particulate 
matter  should  have  no  effect  on  pulsed  propaga¬ 
tion  for  pulse  lengths  less  than  1  nsec  at  both 
1.06y  and  10. 6 p. 
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In  addition  to  these  conclusions  regarding  the  phenome¬ 
nology,  one  can  conclude  that  simultaneous  measurement  of  pulse 
transmission  and  plasma  emission  produces  an  accurate  means 
for  determining  the  breakdown  threshold  of  gaseous  as  well  as 
particle-laden  atmospheres. 
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APPENDIX  A 

DETAILS  OF  TIIE  14ETHOD  OF  QUASILINEARIZATION 
FOR  APPLICATION  TO  THE  LSC  PROBLEM 

A.i  DERIVATION  OF  EQUATIONS 

The  theoretical  basis  underlying  the  method  of  quasi¬ 
linearization  has  been  presented  in  Section  2.4.  The  govern¬ 

ing  differential  equation  to  which  the  method  is  to  be  applied 
is  Equation  (22a) .  For  simplicity,  we  define  X  2  A  +  A  and 

C  IT 

j  ^c+^r*  rrllus'  Equation  (22a)  can  be  written  as 

,T/dT  m\  _  d2T  _  -ldX  /dT\2  ,  PoVo  dh  dT  kL!}>"j  ,, 

Idx'1/  dx2  XdT  \dx/  +  X  dT  dx  X  '  (A-1) 

over  the  domain 

-°°  <  x  <  00  . 


Equation  (A-l)  is  a  second  order  equation  (or  equivalently 
two  first  order  equations);  however,  we  can  apply  the  method 
of  quasi-linearization  (Equation  22)  without  separating 
Equation  (A-l)  into  two  first  order  equations;  viz, 


d2Tk+1 


dlt  /mk+1 

dT  V1 


(A-2 ) 


In  order  to  simplify  the  notation,  we  define 


dll _  _  poVo  dh  o  1  dX_  dT 

d^dTj  X  3f  X  dT  dx 
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Vhen  for  the  (k+l)st  iteration,  Equation  (A- 2)  can  be  written 
as  a  linear  second  order  differential  equation 


d2Tk+1  k  dTk+1 

T~2  P  dx 

dx 


k 


-  r 


Tk+1  +  Sk 


=  0 


(A-5) 


-oo  <  X  <  00 


where 


k 

s  = 


k  d' 


,k 


dx 


k„k 
r  T 


-  II 


k 


(A- 6) 


Equation  (A-5)  should  be  integrated  and  the  boundary  condi¬ 
tions  4,  5  and  6  applied  at  the  far  upstream  (x  =  -°°)  and 
downstream  (x  =  °°)  boundaries.  In  order  to  satisfy  the  boun¬ 
dary  condition  at  both  ends,  we  first  integrate  Equation  (A-5) 
from  x  =  0  to  -»>  and  then  from  0  to  00 ,  where  the  temperature 

(r:.  )  of  the  wave  at  x  =  0  is  arbitrarily  chosen  (in  many 

m 

cases  it  can  be  chosen  as  the  ignition  temperature) .  Then 
the  wave  velocitv  will  be  determined  by  matching  the  slopes 
of  the  function  at  the  matching  point,  x  =  0. 

In  the  upstream  portion  (0  <  x  <  -00)  ,  Equation  (A-l) 
can  be  transformed  as 
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11  = 
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d  T 


dx 


,2 


1 

A  dT 


/  dT_\2  PoV 

\  dx '  /  A 


o  dh  dT 
dT  dx  1 


k  >- j 
c  J 


X  '  <“ 


where  x'  =  -x.  Obviously,  by  comparison,  Equation  (A-l)  can 
be  used  for  both  regions  with  negative  velocity  (— v  )  for  the 
upstream  portion  and  positive  velocity  (v  )  for  the  downstream 


portion .  Then  Equation  (A-5)  is  valid  from  -°° 
the  appropriate  velocity  for  each  region, 


x  <  oo  with 


d2T 

dx2 


"  P 


dT 

dx 


-  rT  +  s  =  0 


( A-7 ) 


0  <_  x  <  00 

where  superscripts  k  and  (k+1)  have  been  omitted  for  conven¬ 
ience  . 

To  apply  the  method  of  quasilinearization,  it  is  re¬ 
quired  that  a  particular  solution  of  Equation  (A-7)  be  inte¬ 
grated  from  x  =  0  to  «  with  the  initial  condition 


T 

P 


0  at  x  =  0 


Then  two  homogeneous  solutions  of  the  homogeneous  equation 


P 


dT 

dx 


rT 


0 


(A- 3) 


are  obtained  with  the  following  two  sets  of  initial  conditions 

at  x  =  0 

at  x  —  0 
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By  virtue  of  the  properties  of  a  linear  equation,  the  total 
solution  should  be 


T  = 


T  +  C,T  +  C-,T, 
p  1  c  2d 


q  =  -r—  =  a  +  C, a  +  C  a 
1  dx  Mp  l^c  2^d 


(A-9 ) 

(A-10) 


I.'ow  we  can  apply  the  two  point  boundary  conditions  to  deter  ¬ 
mine  the  two  constants,  and  C^.  Applying  boundary  condi¬ 
tion  T  =  T.  at  x  =  0  to  Equation  (A-9)  gives  C,  =  T.  ,  and 
in  dT  1  in 

applying  boundary  condition  q  =  ^=  0atx  =  °°  determines 

C9  as 


V">  +  Tin<1c(") 
l2  =  -  - 57F5 - 


Once  dT/dx  =  0  at  x  =  -®  is  satisfied,  the  other  boundary  con¬ 
dition,  T  =  T  at  x  =  -«,  is  automatically  satisfied  by  Equa¬ 
tion  (22).  In  practice,  instead  of  integrating  Equation  (A-5) 

dT 

to  infinity,  one  integrates  to  a  large  x  where  ^  approaches 
zero . 

Since  the  method  of  quasilinearization  is  an  iterative 

•  •  •  •  V 

scheme,  an  initial  trial  solution  is  guessed,  then  p  (x) , 
r'(x)  and  s  (x)  are  calculated  according  to  Equations  24,  25 
and  27.  Due  to  the  drastic  changes  at  the  wave  front,  it  is 
found  that  the  integration  of  Equation  (A-5)  in  x  can  be  best 
performed  by  the  following  simple  method.  First,  we  divide 
the  x  coordinate  from  0  to  x _  into  many  intervals  with  dif- 
fetent  spacing  h.  The  value  of  h  depends  upon  the  changes 
of  p(x),  r(x)  and  s  (x)  with  respect  to  x.  If  h  is  small 
enough  p,  r  and  s  can  be  considered  as  constants  in  each  inter¬ 
val.  Then  Equation  (A-5)  becomes  a  second  order  linear  dif¬ 
ferential  equation  with  constant  coefficients  over  each  inter¬ 
val 
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d2Ti+l  ,  ,  dTi+l 

“d—  -  P<Xi)  -J —  -  r<Xi-'T;+1  +  S(.-:i)  -  0 


(A- 12) 


Tho  solutions  of  the  above  equations  are  straightforward  and 
depend  on  the  values  of  p(x.)  and  r(x.).  The  solutions  are 

aih  a_h  s  (x .  ) 

4-  z  .  1 


rTxTT 


where 


uul 


A 


D 


T  dTi  a_s(x.)~|  / 

=  [-  sr  +  “2Ti  -  ihrr-J  / 


(°t2  a]_)  / 


=  „  als  (xi)  1  / 

.dx  alTi  +  ~r  (x.)  J  / (a2  “  ai 


|Ul|  =  1 

|a2  j  2 


p(*i)  1  yp(x.)2  +  4r(xi) 


(A-13) 


For  P(x.)  +  4r(x..)  >  0;  and, 


m/v  »  _  ah.,  s(x  ) 

iUi+lJ  "  °  (AcosBh  +  BsinBh)  +  — 1. 

r  (x7) 


where 


A  -  Ti  -  s (xi)/r(xi) 


B  = 


1  idTi 

£  j  dx  a 


T.  - 

l 


s  (x  ) 


rTxTT 


l  -> 


a  = 


2  P(xi> 


and 


“  2  ^~\_P  ^xi  ^  2  +  4r  (x^ ) 


(A— 14) 
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for  p(xi)  +  4r  (x^)  <  0.  For  the  homogeneous  equation  (A-8b) , 
the  above  solutions  are  valid  with  s(x^)  0.  After  each  inte- 

(jr  1 1  i  on  with  respect  to  x,  all  T  (x.),  q  (x  ),  T.(x.)f  q  (x .  )  , 

VV  and  values  are  Stored.  Then  the  constants  C1 

and  C ~  are  determined  from  Equation  (A-ll).  Finally  the 
k+1  k+1 

value  of  T  (x^)  and  q  (x^)  ato  calculated  according  to 
Equations  (A-9)  and  (A-10)  and  they  are  ready  to  be  substi¬ 
tuted  into  Equations  (A-4)  (A-5)  and  (A-6)  to  calculate  the 

functions  p,  r  and  s.  The  iterative  procedure  continues 
until  a  convergence  criterion  is  reached  for  the  constant  C2 
in  two  successive  iterations. 


A. 2  METHOD  OF  SEGMENTATION 

The  principal  numerical  difficulty  involved  in  solving  a 
two  point  boundary  value  problem,  especially  where  the  "other" 
boundary  condition  has  to  be  satisfied  at  a  large  value  of 
x«,  <or  at  infinity)  is  numerical  "overflow"  of  the  computer 
during  the  process  of  integration.  For  example,  consider  the 
following  equation: 


R  +  S 


0 


(A-15) 


with  boundary  conditions  specified  as 
y  =  yo  at  x  =  0 

^  =0  at  x  =  00 


(A-16 ) 

(A-17) 


The  quantities  P,  R  and  S  are  all  positive  constants.  Assume 
2 

(P  -  4R)  is  positive,  then  the  solution  of  Equation  (A-15)  is 
u  x  a2x 

y  =  Ae  1  +  Be  +  S/R  (A-18) 
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whet  e 


l°il 


°ii  i  r, 

•2»’  5  L 


-  4R 


Obviously,  is  a  postive  constant,  whereas  a 2  is  negative. 

After  applying  the  boundary  condition  at  x  =  «>,  the  constant 

A  is  set  to  zero,  and  the  exponentially  growing  term  e 

will  not  effect  the  solution.  However,  this  procedure  does 

not  apply  for  a  numerical  integration.  Due  to  the  fact  that 

the  boundary  condition  at  x  =  00  (or  x  =  xb)  will  not  be 

applied  until  after  the  numerical  integration  from  0  to  xro 

is  finished,  obviously,  the  positive  a.  will  cause  the  inte- 

■l*  c XX 

gration  to  be  terminated  whenever  the  term  e  1  is  larger 
than  the  upper  limit  of  the  computer.  Hence  the  integration 
cannot  be  continued  and  the  problem  cannot  be  solved.  This 
is  a  very  well  knowr  problem  in  numerical  integration. 

In  dealing  with  this  problem,  we  propose  that  the  fol¬ 
lowing  procedure  be  applied  to  the  method  of  quasilineariza¬ 
tion.  First,  we  divide  the  entire  integration  path  into  n 

segments  which  are  chosen  such  that  within  each  segment,  the 
on (xi+i“Xi) 

term  e  is  well  under  the  upper  limit  of  the  com¬ 

puter.  Then,  two  homogeneous  solutions  and  one  particular 
solution  can  be  obtained  for  each  segment  with  the  method 
described  above.  The  total  solutions  in  each  segment  are: 


T.  =  T.  +  C0 .  .  T.  +  C0.  T . , 
1  ip  2i-l  ic  2i  id 


i  =  1, 2,  •  •  *n 


qi  “  qip  +  C2i-1  qic+  C2i  qid 


(A-19 ) 
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5  m 

l- 


with  initial  conditions 


•  (X.  ) 

ip  1-1 

=  o  ; 

‘iip(xi-i) 

=  0 

ic<xi-l> 

=  1  ; 

lic(xi-l> 

=  0 

i  =  1 , 2, • • *n 

(A-20) 

id<xi-l> 

=  0  ; 

‘Jid<xi-1) 

=  1 

• 

AC ter  each  segment  has  been  evaluated,  the  two  point  boundary 
conditions  can  be  applied  to  the  systems,  i.e., 

T,  (x  )  =  T,  (0)  =  T. 

1  o  1  in 


q  (x  )  =  q  (x  )  =  0 


(A-21) 


At  the  same  time,  matching  conditions  are  also  applied 
each  joint  of  every  segment; 


W  -  Ti+i<xi> 


i  =  1,2, • • *n-l 


(A-22 ) 


q  .  (x . )  =  q  .  ,  ,  (x . ) 

^i  i  ni+l  i 

Therefore  there  will  be  2n  constants  C-^,  C2'*'('2n  to  ^eter- 
mined  by  the  above  2n  algebraic  equations.  With  some  simple 
manipulation,  Equations  (A-21)  and  (A-22)  can  be  written  as 


C,  =  T. 

1  in 


(A-23) 
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imipiu y  IVFH ' jpn»i.ij!^p  i  .■» 1 » 1 


C2 

U2  A2C1 

a3(b3-i) 

C3 

°3 

a4  <B4-n 

C4 

• 

• 

D4 

• 

• 

• 

• 

A2n-2  ^B2n-2~1^ 

C2n-2 

D-  n 

2n-2 

C2n-1 

B2n-1 

A2n  B2n  . 

_C2n  - 

+C_  , 

L  2n  2n+l_| 

(A-2  4 ) 


where 


At  =  T.  (x.  ) 
2i  icv  i' 


A 


2i+l 


q  .  (x .  ) 

WV  -  Tid(xi>  i^ry 


1C  1  ' 


B2i  "  TidUi> 


i  =  1, 2 ,  •  •  -n 


B  =  gi°!Xl| 

21+1  Tic(Xi> 


u2i  * 


IP  '  1' 


D 


2i+l 


,  x  qic(xi) 

=  -qip(Xi)  +  Tip(x.)  j^y 


1C  1  ' 


Sincu  cquatl°n  (A  24)  is  a  tridiagonal  system,  it  can  be  solved 
simply  and  efficiently  by  the  following  method.  [17] 
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The  basic  formula  in  Equation  (A-24)  can  bo  written 


+  BnC,.,  -  Cra+1  *  ■>»  •  2  1  i  2" 


Nov/,  let  us  assume  a  solution  of  C  has  the  form  as 

m 


(A- 2 5 ) 


C  =  v  C  ,  .  +  g  , 
m  ra  u+1  Jm  ' 


(A-26a) 


or  for  C  , , 
n-1 


Cm-1  -  wm-lCm  +  gm-l  *  (A2Gb) 

Eliminating  Cfn_1  between  Eqs.  (A-26a)  and  (A-26b)  gives 


•,  D  -  A  g  . 

q  -  _ 1 _  r  ,  m  r.vna-l 

u  A  v/  +  B  m+1  A  v/  ,  +  B 
n  m-1  m  m  m-1  n 


Obviously  the  and  g  are 


( A-2  7 ) 


A  W  , 

+ 

B 

n  m-1 

m 

D  -  A 

a 

n  ] 

ivm-i 

A  w  , 

+ 

B 

m  m-1 

m 

(A-28a) 


q  =  — - -  -  A 

Jm  A  w  ,  +  B 
m  m-1  m 

To  start  the  calculation  procedure  one  first  observes  that 
Cf  -  v/l^2  +  gl 

however  since  from  Equation  (A— 23) ,  we  know  that 


C,  =  T. 

1  in 


(A-28b) 
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It  is  apparent  that 

Wi  =  0  and  =  Tin  .  (A-29) 

ihen  all  the  w's  and  g's  can  be  calculated  easily  provided 
that 


A  W 
11 


m- 


_1 

1 


+  D 

m 


<  1 


V/itn  tiie  definition  of  A^  and  B^,  we  can  show  that  the  above 

criterion  is  always  satisfied.  Once  all  w's  and  g's  are 

known,  it  is  simple  to  compute  all  constant  C  's  according 

to  Eq.  (A-26a)  with  C„  ,  =  0. 

2n+l 

To  demonstrate  tne  above  method  we  choose  to  solve  the 
following  nonlinear  second  order  ordinary  differential  equa¬ 
tion 


'2 

dx^ 


“  2y)  -  b  +  ay2 


with  boundary  conditions 


7  =  yQ  at  x  =  0 


(A-3Q) 


iz  = 

dx 


0  ,  at 


x 


00 


where  a  and  b  are  positive  constants.  This  equation  can  be 
solved  exactly by  defining  a  new  dependent  variable  u  a 


u 


Then  Equation  (A-29)  can  be  written  as 


(A- 31) 


du 

dx 


uu 


-b 


(A-32) 
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Integration  of  equation  (A-31)  gives 
u  =  c^  e  +  b/a  , 


or 


dv 

dx 


ax 


(A-33) 


duo  to  the  boundary  condition  at  infinity,  C1  is  zero.  There¬ 
fore  the  exponential  grov;ing  tern  eaX  is  eliminated.  Integra¬ 
tion  once  more  of  Equation  (A-33)  gives 

y  =>/Ftanh  (c  +Vl *) 

After  applying  the  initial  condition  the  final  solution  is 

Y-  JE^yo  +  tanh^X  (A--34 

1  +  -if  Yo  tanh  >/Fx 

We  plot  the  above  solution  for  the  case  where  a  =  10,  b  =  40 
and  yQ  =  1  in  Figure  A-l.  Also  shown  is  the  first  few  itera¬ 
tions  of  the  numerical  solution  by  using  the  method  of  quasi¬ 
linearization.  notice  that  since  a  =  10,  the  exponentially 

lOx 

growing  factor  is  e  .  Therefore  in  the  numerical  calcula¬ 
tion  the  total  integration  path  (from  0.0  to  10.0)  has  to 
be  divided  into  5  to  8  segments  (the  number  of  segments  is 
determined  automatically  by  the  computer  program  whenever  the 
value  of  integration  exceeds  a  prescribed  limit) .  Without 
segmentation,  the  integration  can  only  proceed  to  about 
x  -  1.5  without  an  "overflow".  However,  as  one  can  see  from 
Figure  A-l,  the  slope  of  the  first  integration  will  not 
approach  zero  unless  the  integration  is  carried  beyond  x  =  6.0. 
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Therefore  the  method  of  segmentation  is  very  important  in 
applying  the  method  of  quasilinearization  to  some  nonlinear 
equations.  Also  it  is  shown  in  Figure  A-l  that  in  general 
the  method  of  quasilinearization  is  a  powerful  method.  With 
an  un-inspired  guessing  of  the  initial  trial  function,  after 
only  three  iterations,  the  numerical  solution  coincides  with 
the  exact  solution  to  within  a  half  of  one  percent, 

A .  3  hUTlIhllCAL  1ETII0D 

Because  of  the  highly  nonlinear  coefficients  in  Equa¬ 
tion  (A-2)  (see  Section  2.3),  certain  precautions  must  he 
taken  before  routinely  applying  the  method  of  quasilineariza¬ 
tion.  Since  an  iterative  procedure  is  used,  the  intermediate 
solution  will  either  overshoot  or  undershoot  the  exact  solu¬ 
tion.  It  is  found  that  excessive  overshoot  or  undershoot 
will  cause  a  termination  of  the  numerical  scheme.  With  a 
large  value  of  temperature  (W . ) ,  A  can  cause  "overflow" 
because  of  the  exponential  fitting  of  A,  Equation  (14). 

On  the  other  hand,  a  negative  (although  it  is  physically 
meaningless,  it  can  occur  in  an  intermediate  stage  of  an 
iteration  scheme)  will  stop  the  calculation  duo  to  the  unde¬ 
fined  operation  of  raising  a  negative  number  to  a  partial 
power,  sucii  as  the  calculation  in  Equation  (9).  To  avoid 

these  problems,  the  calculated  temperatures  are  forced  to 

4 

exist  in  the  range  of  10  £  T^/10  £  0.  At  the  same  time, 

special  care  must  be  exorcised  on  the  initial  trial  function. 
Another  problem  which  existed  in  the  iterative  procedure  is 
the  zoning  problem.  It  is  found  that  in  the  iterative  process, 
the  slope  of  the  temperature,  dT^/dx,  can  change  violently 
from  one  iteration  to  another  before  it  approaches  the  final 
solution.  Due  to  this  change,  the  values  of  a  (Equation 
(A-13)  and  (A-14))  will  become  exceedingly  large  such  that 
i  iu'  integration  of  q^  or  is  terminated  due  to  "overflow" 
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within  one  integration  step  (this  can  not  ho  helped  by  the 
method  of  segmentation,  since  a  segment  consists  of  more  than 
one  step) .  Therefore  a  rezone  capability  is  built  in  the 
numerical  program  to  reduce  the  size  of  the  integration  step 
to  continue  the  integration. 

As  we  mentioned  before,  because  of  the  above  numerical 
difficulties,  the  initial  trial  function  must  be  carefully 
selected  in  order  to  assure  a  fast  convergent  solution.  To 
demonstrate  the  method  of  selection  of  the  initial  trial 
function,  we  './ill  describe  the  procedure  to  obtain  a  solution 
Cor  the  case  of  Bunk in s ^  experiment.  A  neodymium  laser  is 
incident  on  air  with  spot  size  .1  -  0.12  5  ci.i  and  the  absorbed 
flux  of  k_  -  1.57*10^  watts/cm3 .  ^  To  obtain  a  solution, 

Ii 

',/e  have  to  construct  the  initial  trial  function  from  a 
simpler  situation  where  the  solution  can  be  easily  obtained. 
Lot  us  rev/rite  the  correlation  function  of  X;  Equation  (14) 

as 


=  2 . 2 * 1 0 8  ea(T/10  1*75)  T  <  17500  °K 


=  2.2*108  e3(T//1°  ~1,75)  ,  T  >  17500  °K 


(A. 35) 


where  a  =  6.4.  The  first  step  consists  of  using  an  arbitrarily 
chosen  exponential  function  as  the  initial  function  for  ten- 

g 

perature  to  solve  the  case  where  kT £  =  1.57*10  ,  a  =  3  and 
neglecting  the  conductive  energy  loss,  i.e.,  j  =  0  (or  R 
equal  to  a  large  number,  for  example  R  =  12.5  in  Equation  (3)) 
with  a  given  velocity  v  .  In  general,  it  takes  about  eight 
to  twelve  iterations  to  obtain  a  solution  with  less  than  one 
percent  of  error  in  two  consecutive  iterations.  Then  by 
using  the  final  solution  of  the  above  case  (a  =  3,  R  =  12.5) 
as  the  initial  trial  function  we  can  obtain  a  solution  within 
two  to  eight  iterations  for  the  case  where  a  =  3,  R  =  5.  By 
the  same  method,  we  then  proceed  to  compute  solutions  for 
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cases  where  a  =  3,  R  =  1.25;  a  =  3,  R  =  5;  a  =  3,  R  =  .25; 
a  =  3,  R  =  .125;  A  =  5,  R  -  .125;  a  =  6,  R=  .125;  and 
finally,  for  Dunkin's  case  vhere  a  =  6.4  and  R  =  .125.  Usually 
Usually  two  to  four  iterations  are  needed  for  each  case,  ex¬ 
cept  when  the  iterated  solutions  oscillate  between  two  Units 
with  a  very  slow  convergence.  Then  r.n  average  method  is  used 
and  usually  it  takes  about  four  to  eight  iterations  to  converge 
to  a  final  solution.  Instead  of  using  the  ith  iterated  solu¬ 
tion  (T.)  as  the  trial  solution  for  (i  +  l)th  iteration,  an 
average  method  using  (T.4T._^)/2  as  the  trial  solution  for 
T.+^  is  utilized  in  order  to  expedite  the  convergence  process, 
bach  iteration  takes  about  four  seconds  of  computation  time 
on  the  UI1IVAC  1108. 

As  we  dr-cribed  in  Section  2.2,  the  LSC  propagation 
problem  is  an  eigenvalue  problem.  The  propagation  velocity 
must  be  determined  by  matching  the  upstream  and  downstream 
;olutions  at  x  =  0.  All  of  the  above  described  numerical 
methods  can  be  applied  to  either  the  upstream  solution 
(negative  v^)  or  to  the  downstream  solution  (positive  v^  ) 
with  a  given  velocity.  Then  by  varying  the  velocity,  the 
gradient  of  the  temperature  at  x  =  0  will  eh  inge  accordingly. 
Finally,  a  graphical  matching  (Figure  A- 2)  (or  a  numerical 
matching  on  the  computer)  of  the  slope  of  the  upstream  and 
downstream  temperature  profile  determines  the  propagation 
velocity.  With  the  velocity  known,  we  then  calculate  the 
matching  temperature  profile  of  a  propagating  LSC  wave. 
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APPENDIX  B 

USER'S  GUIDE  TO  PROGRAM  TEMP 

TEMP  is  programmed  in  FORTRAN  IV  and  run  on  tue  UUIVAC 
1103.  It  is  a  double  precision  program.  This  program  is 
designed  to  solve  the  temperature  profile  of  a  one-dimensional 
steady  state  laser  supported  combustion  wave.  The  numerical 
method  used  in  this  program  is  the  method  of  quasil inear iza- 
l.iou  with  the  method  of  segmentation  (see  Appendix  A)  .  TEMP 
calculates  the  downstream  temperature  profile  as  well  as  the 
upstream  temperature  profile  depending  on  the  sign  of  the 
given  velocity.  For  a  positive  velocity  the  downstream  tem¬ 
perature  profile  will  be  computed  whereas  a  negative  velocity 
will  give  tae  upstream  profile.  In  order  to  obtain  rapid  con¬ 
vergence,  usually  the  initial  function  for  any  iteration  is 
the  final  solution  of  iteration  (see  /Appendix  A)  .  Therefore 
a  tape  is  required  to  store  the  final  solution  for  each  case. 

In  the  following  paragraphs,  wo  first  define  the  input 
variables  and  output  variables.  Then  a  typical  run  with 
input  and  output  will  be  displayed.  Finally  a  listing  of 
TEMP  will  be  included. 

IUPUT  ViV.g/VBLES 

.IEThRT  an  integer  switch. 

0;  the  initial  function  is  arbitrarily  set  in 
the  program 

1;  the  initial  function  is  read  from  the  tape. 
JULOCK  an  integer  number 

When  MSTART  =  1,  the  initial  function  will  be  read 

from  "JBLOCK"  th  block  in  the  tape. 
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YId 

mo 

u 

ah 

TiiUX 

RyvJIUS 

i:ps 

output 

i 

HOLD 

.  V 

TP 

Jp 

TC 

>c 


a  double  precision  nui.iber. 

Die  ignition  temperature  or  a  starting  temperature 
at  x  =  0  in  10  °K,  usually  Till  =  1.2. 

a  double  precision  number. 

T!ie  density  of  air  at  far  upstrean. 

XI 10  =  1  .18  *  10“3  g/cr.v* . 

a  double  precision  nunber. 

The  cjiven  propagation  velocity  (cn/sec) 

U  >  0;  downstream  temperature  profile  is  calculated. 
U  <  0;  upstream  temperature  profile  is  calculated. 

a  double  precision  number. 

exponent  of  correlation  formula  of  X  (Eq. 

a  double  precision  nunber. 

Vac  quantity  ];$  (crg/cm3  sec) 

a  double  precision  number. 

Tiie  radius  of  incident  laser  bean  (cn)  . 

a  double  precision  number. 

Convergent  criterion 

VAXIAbLi.S 

All  variables  are  in  double  precision,  except  IC. 

Constant  C2 

Constant  C2  in  previous  iteration 
X  coordinate  (cm) 

Particular  solution,  T 

p 

Particular  solution,  q 

P 

Homogeneous  solution,  T 

Homogeneous  solution,  q 

*c 

Temperature  (10^  °K) 
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4 

>7  Yoi iperature  gradient  (10  °K/cn) 

P  Function  p(T,  dT/dx) 


Function 

r  (T,  d 

T/dx) 

G 

Function 

s  (T,  d 

T/dx) 

i:i 

Ll  =  cx^ , 

when 

k  =  0 

L!l  =  a , 

when 

k  =  1 

1:2 

b2  =  a2, 

when 

k  =  0 

E2  =  6, 

when 

i — i 

ii 

iC  An  integer  switch 

K  =  0,  exponential  type  solution 
K  =  1,  sinusoidal  type  solution 

Table  A.  1  shows  a  typical  input  data  of  TEMP .  Jotice 
that  the  initial  trial  function  is  read  from  the  9th  bloc!; 
ui  the  tape  (see  Section  6  for  detail  of  numerical  calcula¬ 
tion)  .  Table  A. 2  shows  the  output  of  the  calculation  for  the 
input  data  shown  in  Table  A.l. 


*-2  344 


SIM^UT 

MST4RT 

J8.0CK 

T I  h 

RH3 

U 

•  X 

FLUX 
ROI  US 
FP  S 

*r  no 


TABLE  B.l 


TYPICAL  INPUT  DATA  FOi*  TEMP 


Reproduced  from 
best  available  copy. 


♦  1 
♦9 

.  i?orcoLoooiuijori  cgo*cc  i 
.  ;  130003 30 03 130 13  30  0-002 
.  1  S0CQQ0Q0Ql)0Q0CL'C03+LCt< 
.  53399999999  9  9911990  *001 
.  iOOQOOOOOOOUGGU  LU3*  0  l  ! 
.  *50000300003303030  0*300 
.  lOOOGOGGOOCGLGCl  LCD- OC  l 
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CALL  N  T  R  A  0  I  !  l  »2»|0»J*LU»22* 

IF  (  td  «  N  L  «  I  Q  )  *  E  T  U  H  N  -j 
P‘4  CO'ITINUC 
23  CUNT i N  Jf 
P  »  I  N  T  l  ,  1  t 

iCI-M  (  ‘O')  To*  '  ’  »  /  (  I  H  ,  1IM0LI-.4)  ) 
f  0  n’  M  m  T  (  *  C  T  =  ’  »  /  (  1  H  t  IP  10L'IU«4  )  ) 

print  13 

XL  u  I  n‘1 1  'J'o  of  |  LOUP  uN  T  £  Ml’  f  h  a  1 ,J*  F 

■>K  1  'T  3U,  c  * 

3-3  F  (;  H  M  M  T  (  I  H  I  ,  •  C  X  “  ’ 

!  C  r  a  . 

cu  ,1  i  •;  o  (■' 

u  t« i c  r ♦  i 
i pf  * i pp 
ki  » i 

ULI1)«I 

no  2 j  i * i »  ip 

IF  (  r  E  <  I  )  .LL.  3.  .OH  .  TE  I  1  )  •  Cl  T  •  3  •  D  t  )  (,0  TO  86 

NO  T  U  8  7 

-3  8  TE(I)*TC(1-|) 

4  V  C  U  f  J  T  I  N  U  r 

20  CO  IT  1 NUE  ' 

T‘,«.  ^3 

cu-w'  jtf  l>t»?  i  . /.  r  i  ;f.s  d*ot  ,02a0T  .uaoT  ,p?bot  tncnT  ,  ano  p 

JO  S-  I  a  1  ,  I  f 
IT  *  I 

call  0  L  » 1  (  n  (  I  >  ,  0  T  0  A (  i  )  ) 

p  =  1  .  /  I 
om2»u4«*2 

!’(  :  » «  o  n  •  t  i « t  -4«0*oaot-2.*,jbi)T*otux(  j  >  > 

w  «  I  I  *  0  T  0  A  I  I  )*oU2*(c,*(B*02MDT-DAoT*DBnT)«l.E:-3-(b*D2RUT-Otj!>T**2>*ijT 
1  0  *  I  1  I  )  -  |  I)  C  0  T  •  0 .  j  -  L  •  0  B  2  •  0  8  0  T  )  •  l  •  t  -  H 

<  I  )  *P  I  I  >  •  0  T  U  A  (  J  >  (  j  )*TC(  1  )  —  0  tJ  •  C  l»F>***G*L)AjT-lJ,iUT*UTr)X{  I  )  )  *OTOa  (  I 

1  •♦C/wM.L-4 

13  format  (in?,  *xx,  tl,  jaDT,  la'Tdac,  olcot,  l>2lloT,  lamdaw,  dlrdt, 
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in?L«uT  I,  fl»  ObDT,  i)  2  j  T  ’  /  i  X  i  ’  I)  T  u  X  ,  FXL»  bKL»  0 1*'  L  0  T  »  PHl,  ’> 
l  OJOT,  AJC,  UJCOT,  C,  JCOT*//) 

!OifAr  <  i  2  <  •  /■  .  ii:.m)/2x,  1 1  (diu.h.  i  x  i ,  n9.->> 

COnTi NOE 


!  r.  TLgRATl  TO  OJTAlU  DTDX 


*  >  l  1  i  »  J  . 

JP ( 1 ) «0 • 

T  P l I  I • 
r4  C  l  1  )  ■  0 . 

TC (  I  )  ■  1  . 

OD (  1  i • I • 

TO ( 1  I ■ 

N  1  a  1 

QO  h  l  M  a  2  ,  J  N 
N  ■  IJ  N  l  H  )  -  l 
•1»MH  (  M  ) 

)  0  6  „  1  c  N  1  , 

I  1  3  1  ♦  I 

a  »  X  «■  ,i 

3w«St  1  )/K(  l  > 

S  a  L  *  (’  (  1  )  •  •?*4  •  (  I  » 

H  IjAU  6  3,63.0** 

HF  T A«.b*bOKT  <-S\U ) 

\  A  ■  •  G*P (  1  ) 

UMl'l 

I I  l  l  l  )  *  A  A 

1  E  2  {  i  )  ■  H  F.  T  A 
*IK  L-HKL*  1 
K  A  «*  E  a  P  (  A  A  *  M  I 
OA«Bt.TA*H 
CHAsCOS  (  ti  A  ) 

S  e  A  «  b  I  N  (  H  A  » 

C  1  *  T  L  (  I  ) 

C2«(0CtI)-AA*CI)/BETA 

I)  1  «TU  (  I  ) 

0  2«(uDU)-AA»L’U/BETa 

P 1 »TP ( 1 ) -SW 

i»  2  3  (  N  P  <  1  >  -  A  A  •  P  t  >  /  B  E  T  A 

NF’  (  I  1  )  *  E  A  ♦  (  (AA*P|^«ETA«P23*Ct3A+(AA*p?-i»ETA*IM  )  •  b  b  A  ) 
IF  ( A  b  S ( 0  P (  1  1  )  I  . G  T .  l.OU  )  GO  TO  74 
TP  (  J  1  )  «r  A*  (  P  1  •CnA*P2*S|j  A  ) 

•3  Cl  I  j  I  *  E  A  •  <  (AA*Ci*BETA*C2<*tilA*(AA*C2-bETA*Ci  )*SbA» 
TCI  1 i  )«EA*(Cl*CBA*C2*SbA) 

D  <  I  l  )  »  F  A  •  <  (aA*D1*b£TA*D2)  •CliA*(AA*02-tjETA«01  )  •  S  B  A  J 
T  lj  (  I  l  )  •  C  A  *  <  *>  1  •  C  ri  A  ♦  D  2  •  5  B  A  ) 

F  0  I?  M  A  T  (IX,  6  (  ]  X  ,  D  I  I  •  *♦  )  »  I3X,  U  I  1  t  4  ) 

0  TO  6tJ 
?  A  L  »  j  0  P  T  (Sal) 

a  l I • . 5  * ( P  <  I  ) ♦SAL  > 

aL2**5*<P<1>-5al> 

1 i  (  I  )-0 

K  i  (  i  )  *  A  0  t 


il-2344 


It  2  (  i  >  *  A l  7 

I 

A  ?  I  *  A  l  2  -  A  1  I 
A  L  I  H  ■  A  L  I  •  m 
t  A1*LXP»/.H»H) 
t  A?»LaF(AL7«n) 

/L>  C  1  51  (  “  U  C  (  I  t  ♦  A  l  2  *  T  !_  I  1  I  t  /  h  2  * 

C2»<UC<I  » -AL  I  • I c  < 1  >  > /A2l 
J  I  *  (  “  U  U  (  I  l  ♦  A  L  7  •  T  I"  (  1  )  )  /  A  2  1 
0  2  ■  (  W  0  (  1  )  •  A  L  1  •  T  o  (  1  I  )  /  A  2  I 
P  1  ■  (  •  0  P  (  1  )  *  A  L  2  •  T  F  (  1  )“AL2*jK)/A21 
P2»(mP(  l  ) “ A |  1 • 7  p (  |  ) ♦AL l *5K ) / A2 1 
JP(I4»»AL1*P|*E;a14AL2*P2*uA2 
IF  {  M h b»  (  u P  (  i  1  )  )  ,gT.  1#dJ  )  GO  TO  7S 
I  / 1>  TP  <  I  I  >  »P  1  A  1  -»P2*E  A2*r>P 

3  C  (  1  i  )  a  A  L  1  *C  1  fCA  I  ♦AL2*C2<LA2 

r  c ( 1 1 ) *r i *r a »  *c2*ea7 
.’0  l  T  I  isM|M|'tH|*ALi*l)2*LA2 
T  i)  (  1  i  )  «ii  l  •  \  A  1  ♦  u  2  •  E  A  2 
*2  F  n  f*  m  «  T  l  1  A  ,  7  (  I  X  ,  I) !  I  .H  >  ) 

*0  CONTINUE 

ci  1  N  I  «  N 

't1'  I  M  <  KL*  I  )  *  I  (,c 

TC  I  (  r*L  )  =  T  C  (  I  Pt y 
i  C  I  (  N  L  )  *  •<  C  (  I  r  i  ) 

fin  (  m.  )  *  td  » i  pi  i 
ud  1 1  kl  >  »'j:n  i  pi  ) 

TPJ  ( KL  )  .tp (  i Pt  * 

'3  p  1  (  k  L  >  ■  In  p  f  I  p  t  ) 

KL2»KL*2 
CON ( 1 ) »T 1 H 
C  0  N  <  k  L  2  ♦  1  »=;. 

u  COMPUTE  nsr  VALUE  Ks-UP/UC  a  no  test 

73  CONTINUE 

C  0  L  0  *  C  o  N  (  ?  i 

It  (  N  L  •  l  J  <  I  )  go  TO  8  J 

n  i  )*o. 

A G  (  |  )  *COM  I  ) 

XL  1 *KL-  l 
•)0  2 ~  l  i  =  !  .  v  L  l 
1  *  2  -  L 
U  2  =  L  1  ♦  I 
l.?*L  1  -  l 

A  b ( L  2 ) *  0 C !  «L)/rci  ( L  i 
AMLll*TUt  U  ) 

A (LI  J  »  T  C l  (L  t 

A  <  L  2  »  ■  N  O  1  (  1.  »  -  T  )1<L>*AP(l23 
r>  (L  l  >  *-TP  1  <L  > 

0(L2)»-0PI  t  L  » ♦ r  p  1  ( L ) • A  H  (  L2  ) 

'(LI  )■!./(  a(L1  )  •  (LSltAliUl  )  ) 

•''  (  l 2  t  *  l  •  /  £  A  (  l 2  J  (  l  1  )  ♦AM  <L?  >  > 

A  G  (  L  1  )  c  (  U  l  L  I  f  •  A  (  L  l  f  •  A  &  l  L  j  f  J  *  v‘*  (  L  l  ) 
f  i  ■*''»  a  1  >  =  ( (;  (  L2  »  -  A  (  L2  >  •  AG  (  L  J  3  )  •  «  (  L2  » 

1  I  K  l.  2  )  =  '  .  /  <  U  I  i  K  L  l  •  l  K  L  2  ■  i  )  ♦  >3  0  1  (  K  L  i  ’ 
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iT-2  34  4 


*  (i  1  M.  2  )  *  <  -  H  1  t  K  L  >  -'JC  I  lKL)*AG(KLi-l  <  )  *  *  <  X  L  ?  I 
00  2J2  L  *  K  L  2  i  2  •  -  1 

2  22  CON(  l  )  ■«(  l  )  *C0u(  t*l  >  ♦ah  »  L.  > 

GO  T  u  !  8  b 

^1  pKjNT  Ij,  lfGMjm,JK.|,i:L2) 

II  ToRMaT  (  2  (  H  i  l  t>u  l  p  .  4  )  j 

PRINT  6  (.'111  p,  ,  COLO  ,  f*h  I 

FORMATC  <  *  »  .11*012. -t,*  KOLO  -  *  ,  I  P  n  I  2  .  4  ,  lA  .  *  P  p  I  -  •  ,  D  1  2  .  4  ) 

0  0  7  L  J  * !  i K l 
,-1  f  h  ■  1  K  LI  J  ) 

*’  *■  ■  I  K  L  I  J  ♦  I  )-l 
1 2  =  2  •  J 


7 ; 
7C 

Jfc 


fl  6 


*- 


"2 


i  <  I 


I «J2-  1 

>0  7  1  I  *  '1 K.  K  ,  '• , 

Tt(  1  *  *  r  t-*  4  I  >  ♦  C  0  II  j1  I  •  T  c  •  I  )+COn(J2)*T|j(  I  ) 

0  T  •)  X  t  I  j  *  w  P  I  I  )  4  c  0 fi  (  J  ;  )  •  1C  (  I  )  ♦  C  ON  (  j 2  )  •  k l>  (  I  1 
C  0  N  T  i  f,  U  r 
C  0  N  T  j  N  0  £ 

IF  (*Hu  .LO.  lPp)  00  ro  82 
I  P  I  “  I  P  L  ♦  i 

0  0  0  0  1  *  I P  1  I  I  P  P 

T  c  (  I  )  >  r  K  (  I  P  l  ) 


DTDXi I  IsDTOXt  JHQJ 

PRINT  AX(D,  TNII.&lMp.Tcm.QcMI.TEin.oToA,,.  pij,  l(( 
1 1 .  iui  1 1 ,  ir.2(i » ,  u  ( i » 
test  on  * 

IF  UltSl  1C  r/s-icT/H.)  •  G  T  •  l.E-ftl  GO  To  In? 

PUNT  j 

n  0  |  ‘I  i  i*|,  ip,  j  ” 

r/1N!  M  Xa'1''  t,(  1  )  »)P(I  )  iTC(  I  )  *v»C(  I  )  ,Th(  I  ■/  ,|)T(ja  ,  J  ,  ,pj  ,  |  . 

i  ' .  I  f :  i  ( i ) ,  i  r  2  ( j » ,  j  t  <  i  >  ’  ‘ 

CO  IT  j  hoc 


I  4  2 

I  3S 

I  21 
I  78 
I 

1  ‘<3 
I  3  3 
1  OC 


CONI  I NOf 

IF  <ABS<<COII<  2  )-COLD)/COLO)  .GT,  EPS  ,  g0  |0  10 
p  R I N T  1  2  1 

TORMmT  (  3  A  ,  i  j  t  £  0  N  V  £  R  r,  c  S  »  ) 

P  R  I  N  T  3 


00  I  *13  I  *  1  ,  IP,  13 

PnjkT  4  X  X  C  j  )  ,  TP  <T)iOP<II|TQ<iJ#Qc<i),TE(j) 

1  *  i  I  F.  i  <  l  )  ,  I  L  2  <  I  )  ,  i  U  I  ) 

CAtL  nTpan  1  I  I  •  1  I  I Q  .  N , ! Q , 2 2 ) 

IF  1  LO  .  NL  .  I  < j  )  '  LTURn(. 

GO  T  u  2  1 


.DToA<i),p(i),k<i 


Vi  C°N  <  2  )  =  <  -OP  <  I  P  I  -C0f»(  1  )  *Nc  i  I  p  )  )  /  ^0  U  P  j 

I  p  S  K  C ■ K  L ♦ I 

CF  <  KC 1 *  C  0 1 J  »  2 1 

' 1  1  i  A  V  ,  f  0  ,  i  )  r,0  TO  1  7  7 

IF  (NC  .LT,  jj  ]vi  |  0 

EC  <  =f  C~  I 

K  C  2  *  t  C.  -  2 

FC3*nC-3 

IF  <  C  K  <  K  c  >  ,  Ci  T  .  C  F  <  K  C  1  )  )  GO  JO  0  1 
IF  <CK<KCI)  .LT.  C  K  (  K  r  2  )  )  V1 0  10  ol 


I  .  S  <  1 


)  .  S  (  I 


>•5(1 


R-2344 


If  IUIKUI  .  (» r .  cmkOI),,,  Tu  rtJ 

CKCK«(CK1K(.)*(X(*CI J )  /  2  , 

if  ( i.KO  *i  j  .  j.\jZ  )  (jV  To  til 

(  A  J  *  1 

I  f  7  '  0  U  T  i  Nuf 

ph  i  r.  I  j  32 

]i?  r0HM-T  Mh'*  c'*.  ’aVlRa&c  comv£kgeot  Mn..nP  is  ub. 
f-CC^CC  ♦  I 

je>  M  ^  I r-  r  i],  ( c  on  c  .jn )  ,  jk*i»  *L2> 
fSlNf  6,  COM  21, COLO 
”0  3  J  J  a  I  ,  K l 
,-1  >  <  ■  .  K  (.  <  J  ) 

**  ■  |  N  L  (  J  ♦  l  J-l 

J2«2 • j 

J  I  *  J  2  -  | 

DO  3  s  I«.NKN,hr 

Tfll»«<TH<|)*(),(j|)*Tc<!,*cuN«J2>*T;>(l)*TL(in/2. 

TTDX  ,  I  >«  (  4H  I  I  ,♦000(01  )  •  4  C  I  1*C0n(J2)»U0(  I  )01[)a  (  I  '  J  /  j 
3  *«  CO  rj  Tin  Up  * 

3  3  coot  ,  no r 

*  ** <J  1  r  I  r  CC  ,  LL  .  f*  I  i,0  TO  3fa 
(  I  7  v 

l/v  fOI^T  (  3>  ,  •  .-F  Koi  MOtVE  aVC«a<.{  i  T  r  •  *  a  T  1  0  r ,  •  j 

o  0  T  :j  !  7  << 

^  If  I  f'  L  •  l  '  i  •  l>  I  M  ,  (i  (J  I  (J  ()  ,4 

II  <  M-  •  L  C.  •  1  i  (,0  To  2  V 

M  K  L  ■  C 

T  f5  I  (  N  L  )  *  (  P  (  I  , 

U  I  ( M  )  «TC  I  |  ) 

I  D  I  (  NL  )  *  I  CM  I  ) 

'■Pi  (  N  L  )  *  C  P  (  I  , 

«  C  I  (  N  L  J  *  0  c  (  1  J 

1  0  I  I  N  L  1  =  U  )  (  I  I 

TP ( I ) «C. 

WP ( I > «C, 

T  C  (  I  ,  ■  1 « 

<JC  I  I  I  «c . 

T  0  I  I  <  «  0  , 

J  f>  (  I  1  *  1  • 

M_«KL*  I 


I  M.  (  K  L  I  a  | 

IP  1  I  E  <  I  I  )  7  b  ,  /  s  ,  7  «s 
•'  ;  of  *  1  - 1 


7 


j 


u 


A  OP, 


bO  To  tis 

°K  I  '•  M  '  1  *  .  Tf’  (  I  )  ,QtM  J  I  ,  (c  (  1  ,  ,oc  , 
1  1  •  I  c  I  I  j  •  1 1 . 7  »  I  )  ,  1 1  (  I  ) 

U  Tui«u  r 


l  ( 


.Oh  M  U  ,P(  !  l  ,i{(  ;  ,  ,s  .  ; 


II".,6>,»1',  I  C  X  ,  •  TP  •  ,  V*,  •  C,  r  •  ,  '//,,• 
l>,'p.»,9.,,P»iyXf.H.|  7X,»S'OX,*LI»,/X.M2‘ 

1  0  r/  ’  ‘  T  '  t-  ‘  3|,  •<(  1  1  X  >  »  2<  E9.S  .  1  X  >  .  I  l  ) 

Fno 

IS  SiMP.SI^P 

SUHHOuTlur.  SUP  (  Y  ,  H  ,  I  ,,  ,  2  i 

IMPLICIT  DOM, jL£  PKtCISim;  {A-H,  0-7, 


r  v- ' , 

,  S  A  , 


9  A  ,  •  kC 

K  *  //  ) 


1  r  x ,  m  • ,  1 0 
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»i  MENS  I  Oi,  Y  (  2V  1  ) 

2  »0  ■ 

1 f  <  »  W)  ?,  3,  M 

o  i  [O,  < .  ,  ;> 

aZ*iY(  I  -  y  i-|  i+n  J  )  )*V 

w  J»  i  ii 
r>  *  A  I  1 1  /  2  « 

M  *  IN/? 
o  =  ij-n 

'  F  l  i;  •  0  T  .  .  I  )  „  J  TO  ; 

3  /  *  2  ♦  i  r  (  i  \  2 

2  wf-  Tum, 


t  \  j 

'  ^  ,  1  o  j  f  t<  1  I  D  L  H  I 

SUHRUU  r  Jilt  0L‘M(U.0Tf)A) 

I  ''’LICIT  OTJdl  ii  PRECISION  <A-t!,U-2> 

OU'JQUt  PRECISION  LaMuaC,  UA'-'UaK,  K  2  1 

''a^a  1  e  r r ^  i  pa  i  2o  i 

CuMiiO,.  nA  )»,Ul'Or,D2AoT,D2j0T,[>CDT,h,C 

inMM'.’oic  r"1'  ‘J"'  oL“noc'!oT,o^cD^2LRoT^,Uu,lI;T,DPu,,AJ,,JCi 

c  0  >1  M  U  u  /  A  A  *  /  A  * 

I  a  T  K 
‘  *  1  •  /  T 


i 


C 

r. 

c 

c 

c 

c 

c 


2  I 


22 


**  •  3  f  31  i  '  •  (  T  /  j  .L4l»*2CP**2/Stc**2 
A  JT=  I  •  0v*22ot.l  2*1 

?  A  0  T  »  j  .  .  1  ?hl  i  i 

M'?mhT(.  r  ad  T  a  •  .  1  f’  E  1  ^  .  'I  ,  *  u  2  A  L>  T  a 


’ » t  Pt i: .h) 


d  *  L  A  h  0  A  C  ♦  L  A  P  J  A  t\ 

LAMDmC-.  I29E7/CJSH(5..(T/1.F.H-|.hs>  >  EPG/SE  OC.l-uEDK 

La  .  !6H7*I  T  /  1  •E'U«*j/n2*LAP(-,J1  )  •  C  u  I  •  •  M  ♦  «* .  .  u  I  •  •  i*  U  .  .  „  .  .  .  . 

♦2*<  •  *0  I  +  2M  •  »  ERr,/SEC-C'1-0LGK  * 

U  I  »  |  2 . 6  5  /  <  T  /  l  ,  E  4  ) 

<2"I*77a(T/|»F  ‘I  >••(-.  7  1  >  c  M  •  •  (  -  j  ) 


IP  IT  .  r,T  .  )  .  7SU3I  DO  Tq  jj 
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APPENDIX  C 

APPROX  I' "1ATE  TREATMENT  OF  RADIAL  EXPANSION 

f  21 1 

According  to  Baker  the  radial  position  of  a  cylin¬ 

drical  blast  wave  is  given  by  the  formula 


Rs  '  <C'1: 

where  R  is  the  characteristic  radius,  ao  is  the  speed  of 
sound  in  the  ambient  air,  and  t  is  the  time,  extrapolated  to 
the  time  that  the  shock  was  hypothetically  initiated  at  a 
line  source. 

If  we  assume  that  the  radial  velocity  of  a  mass  element 
behind  the  shock  is  proportional  to  its  radial  position, 


then 


u(r)  =  ~  .  (C— 3 ) 

The  change  in  mass  in  a  cylindrical  volume  of  radius 
R  and  length  l  is  given  by 

m  =  pu(R)S  /  (C-4) 

where  p  is  the  density  (assumed  constant  in  the  cylinder)  and 
S  is  the  curved  surface  of  the  cylinder.  From  this  formula 
wo  get 


n  =  p  2ttR£ 
=  m/t  , 


(C-5) 
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where  m  is  the  mass  in  the  volume  of  interest. 

The  time  for  a  blast  wave  to  reach  R  after  being  ini¬ 
tiated  at  a  line  source  is 


R 


fco  R  a 
o  o 


(C-6 ) 


The  characteristic  radius  is  given  by  Baker  as 


HE o 

VBYl'o 


R  = 

o  WBYu 


(C-7 ) 


where  Eq  is  the  internal  energy  per  unit  length  in  the  blast 
source,  pQ  is  the  ambient  pressure,  y  is  the  ratio  of  speci¬ 
fic  heats,  and  B  is  a  function  of  y.  For  air,  y  1  5/2,  and, 
therefore,  from  Baker,  B  -  4.  For  a  perfect  gas 


Eq  =  (y-1)  1  kT  R2pq 


/  , . -1  2 
=  (Y-1)  P  TTR  P, 


(C-8 ) 


where  p  is  the  pressure,  T  is  the  temperature,  and  pQ  is  the 

density  of  the  gas  behind  the  shock,  at  the  time  t  ,  i.e., 

o 

when  the  shock  is  assumed  to  be  started  at  R.  Therefore, 


TTp 

R  =  R  f—r  °  .  E- 
o  tly  (y-1)  pQ 


(C-9 ) 


and  from  Equation  (C-6), 


t  =  t  +  t' 
o 


J 


I  i  y-D  12  +  t* 


R 


ao  \  P 


(C-10) 
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where  t'  is  the  is  the  time  after  shod;  is  initiated  at  radius 


The  difference  equations  presented  in  Section  3.4  are 
modified  in  the  following  way; 


1. 

2. 

3. 


The 

has 


continuity  equation  (Eq.  3,10 
m j  added  to  the  right  side. 


of  Chapman) 


The  momentum  conservation  equation  (Eq.  3.1  of 
Chapman)  has  mj+^uj+^  added  to  its  right  side. 

The  energy  conservation  equation  (Eq.  31)  has 
•n+1  n+1  ,  ,  -« 

nm  e^  added  to  its  right  side. 


These  equations  are  then  solved  iteratively  until  convergence 
at  time  tn+^  is  achieved. 

Since  m  is  a  weak  function  of  the  pressure,  it  is 
assumed  that  only  one  shock  is  propagating  from  an  element 
at  any  one  time.  The  shock  is  initiated  whenever  the  ratio 
P/PQ  is  larger  than  some  value  nQn,  and  is  turned  off  when 
the  pressure  ratio  returns  to  some  other  value  noff,  where 
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APPENDIX  D 

A  QUALITATIVE  DESCRIPTION  OF 
LASER- INDUCED  BREAKDOWN 

In  order  to  obtain  a  qualitative  theoretical  descrip¬ 
tion  of  the  breakdown  process,  it  is  necessary  to  consider 
the  equations  which  govern  the  electron  density  and  electron 
energy  gain  and  loss  processes  in  the  laser  field.  It  can 
bo  shown  that  the  electron  continuity  and  energy  conservation 
equations  can  be  written  respectively  as: 


and 


dN 

_ t 

dt 


(v.  -  v  )N  -  DV2N  -  anN2 
1  a  e  e  Re 


dt 

dt 


2  2 
e  E 

.  2  ,  27 

m(w  + 


2ra 

M 


V  -  I  V  . 

c  pi 


Here 

v.  =  ionization  rate 
i 

v  =  attachment  rate 
a 

D  =  diffusion  coefficient 

uR  =  recombination  rate 

E  =  laser  electric  field  strength 

m  =  electron  mass 

gj  =  laser  frequency 

\>c  =  electron  collision  frequency 

I  =  ionization  potential 


(D— 1 ) 


(D-2) 
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N  =  electron  number  density 

t.  =  electron  energy 

In  Equation  (D-l)  ,  the  electron  density  increases  due  to  loniza 
tionsf  and  decreases  due  to  attachment,  diffusion  and  recombin¬ 
ation  processes.  As  shown  in  Equation  (D-2) ,  we  consider  the 
electron  energy  to  increase  due  to  inverse  bremmstrahlung  and 
to  decrease  by  elastic  collisions  with  atoms  (and  ions)  and 
by  inelastic  processes  which  are  here  limited  to  ionization. 

At  the  start  of  the  process,  the  number  of  electrons  in 
the  cascade  increases  exponentially,  i.e.,  the  rate  of  change 
of  electron  density  is  approximately  constant.  Neglecting 
recombination  during  this  early  phase,  it  can  be  shown  that 
Equation  (D-l)  reduces  to 

N  (t) 

~fi -  =  exp{v.  -  va  -  D/A2 } t  =  2k  ,  0-3) 


where  k  is  the  number  of  electron  generations  produced  during 
tine  t.  It  is  common  to  assume  that  breakdown  occurs  once 
some  critical  number  of  electrons,  Ncj_,  is  formed  or  equiva¬ 
lently  after  some  critical  number  of  electron  multiplications, 
kcr'  1S  attained.  Accordingly,  we  find  from  Equation  (D-3) , 


v . 

1 


V  + 
a 


( Jin  2 )  k 


cr 


0-4) 


wuere  we  have  defined  as  the  laser  pulse  tine  required  to 

create  breakdown,  i.e.,  to  generate  II  =  II 

e  cr 

The  breakdown  threshold  can  now  be  defined  by  use  of 
tno  energy  balance  equation  and  requiring  that  during  the 
Lu.>oi  pulse  time,  t  ,  over  which  the  electron  number  density 
has  increased  to  Ncr,  that  sufficient  energy  has  been  added 


107 


R-2344 


to  the  electron  gas  to  create  that  number  of  electrons  and  to 
balance  any  other  energy  losses  which  nay  have  occurred.  Ac¬ 
cordingly,  combining  Equations  (D-2)  and  (J-4),  we  find  for 
the  breakdown  laser,  flux,  F  , 


ncl 


(uSv2) 

c 


lkcr  ln2 


4ircb 


+  V  + 

a 


c 


A‘ 


2nev 

_ c 

i  ;i 
p 


(D-5) 


Here,  we  have  used  the  relationship  between  the  laser  flux  F 
and  the  electric  field  strength  E; 


E  = 


0-6) 


In  order  to  obtain  Equation  (D-5)  we  have  assumed  a  "square” 

laser  pulse  and  an  average  constant  value  for  both  v  and  ~ 

c 

over  that  pulse. 


For  gas  breakdown,  one  can  obtain  limiting  behavior  and 
approximate  scali  ig  information  from  Equation  (D-5).  For 
instance,  for  very  short  laser  pulses,  the  first  term  in  the 
brackets  dominates  the  equation  and  we  find 


F 


I  (w2+v2) 

-R _ si 

v 

c 


(d-7) 


In  this  limit,  breakdown  is  limited  by  the  amount  of  energy 
necessary  to  ionize  the  gas  (kcrI  )  •  The  breakdown  flux  in¬ 
creases  as  ip  decreases.  This  is  the  constant  energy  regime. 

For  long  pulses,  the  energy  loss  terns  become  dominant 
and  v/e  find  that 

(u2+v2) 

r  _ _ c'  e 

£  v  T.  (‘-)~ 

c  loss 
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where  i,  is  some  characteristic  loss  time.  This  is  Lae 

loss 

limit  commonly  referred  to  as  the  constant  flux  limit  of 
breakdown.  In  this  limit,  tiie  flux  to  break  the  gas  down  is 
a  constant  since  the  energy  loss  rate  due  to  diffusion, 
attachment  and  inelastic  processes  increase  as  the  pulse 
length  increases. 

Vii is  simple  picture  of  breakdown  works  fairly  well,  at 
least  qualitatively  for  the  case  of  clean  air.  However,  due 
to  particle  interaction,  the  picture  becomes  more  complicated 
This  is  primarily  due  to  the  fact  that  one  must  now  consider 
the  electron  cascade  evolution  in  a  heated  material,  such  as 
vapor  or  air,  subject  to  rapid  hydrodynamic  expansion  rates. 

In  this  case,  the  energy  loss  processes  themselves  are  a  func¬ 
tion  of  the  pulse  time,  thus  negating  the  long  pulse,  constant 
flux  limit.  Furthermore,  it  is  expected  (and  observed)  that 
additional  losses  depending  on  the  particle  size  and  number 
density  distribution  incorporate  additional  losses  which  result 
in  a  regime  intermediate  between  tne  constant  energy  and  con¬ 
stant  flux  regime. 

The  only  real  conclusion  to  be  reached  at  this  point 
is  that  the  concept  of  simple  scaling  with  wavelength  in  the 
case  of  breakdown  in  particle  laden  atmospheres  is  very  compli¬ 
cated.  However,  there  are  no  experimental  data  available  to 
guide  or  confirm  the  model  at  wavelengths  other  than  10. 6p.* 
Clearly,  such  data  would  be  extremely  useful  in  obtaining 
and  validating  a  model  from  which  simple  scaling  information 
nay  be  derived.  Experimental  data  at  1.06  m  are  presented  in 
Chapter  4. 


Since  the  performance  of  the  work  described  in  Chapter  4,  some 
data  obtained  by  Lencioni  at  Lincoln  Laboratories  have  become 
available  at  l.OGy;  c.f.,  Appl.  Phys.  Lett.  2_5,  15  (1974). 
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